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INTRODUCTION 


The backbone of ballistics has been the spin stabilized projectile. 
Virtually, all ordnance from small arms to artillery has almost 
exclusively utilized the spin stabilized projectile over the last century. 

* 
Its predecessor was the cannonball) and spherical shot. Just as the 
elongated spin stabilized projectile yielded a marked improvement over 
the less efficient cannonball, so also fin stabilized ammunition offers 
great aeroballistic improvements over the spin stabilized projectile. It 
has only been in recent years that the fin stabilized projectile has come 
under serious consideration. Some success was achieved by the Germans 
during World War II with Naval projectile artillery, During the Korean 
War fin stabilized anti-tank ammunition was introduced which improved 
the effectiveness of the shape change because of its low spin, In recent 
years the accuracy of fin stabilized projectiles has improved due to the 
application of the Tricyclic Theory, the use of dynamic supersonic wind 
tunnel tests, and improved launching techniques. Because of their small 


size and the desire for very inexpensive manufacture, the flechette has 


not received the careful attentic-n that it requires to achieve high accuracy. 


It is essential that manufacturing tec] niques, saboting techniques, launch- 
ing techniques, blast suppression techniques, optimized aeroballistic 

design procedures, dynamic wind tunnel tests, accuracy theory studies, — 
computer analysis, and precision firings all be undertaken and optimized 


to achieve good flechette accuracy and low dispersion. 


are 
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The purpose of this study is to explore flechette design and 
performance with a view towards achieving high accuracy and low 
dispersion. Specifically, exploratory firing programs were carried out 
by Frankford Arsenal, by the Ballistics Laboratories and by the 
University of Notre Dame. The results of the Notre Dame Flechette 
Firing Program are summarized in Appendix A. 

A dynamic wind tunnel testing program was also carried out by 
the university on various flechette designs so as to determine the 
essential static and dynamic aeroballistic stability coefficients. The 
results of this dynamic wind tunnel program are summarized in Appendix B. 

f particular importance is the development of a computer theory 
for flechette flight performance, accuracy and dispersion, This theory 
together with an extensive computer analysis is given in Appendix C,. 
Finally, tlechette firi-¢s were carried out in the precision range at 
Frankford Arsenal! and a correlation of theory and experiment is also 
provided in Appendix C. along with a physical evaluation of dispersion. 

Based on the theory, sabot design and lar cher changes were made 


if in order to reduce the values of those parameters which affect dispersion. 


| A second series of firings were conducted and the analysis of the results 


| 


is provided in Apnendix D. 


DISCUSSION 


Exploratory Flechett2 Firing Program 

The exploratory flechette firing program both at Frankford 
Arsenal and at the university have provided an opportunity to measure 
flechette spin, to measure flechette accuracy and dispersion, to 
identify fin damage and body damage due to stripper, to provide an 
approximate measure of dynamic stability at long range, to provide 
a first hand appreciation of the strong blast region and to concentrate 
on sabot design, separation, and transition all as affecting flechette 
flight performance and accuracy. 

In addition a transition ballistic range was set up and optimized 
at Frankford Arsenal to obtain initial condition data using flash x- 
ray photography. A complete description of the set-up is provided in 
Appendix C. 

BRL conducted free flight tests in their transonic Spark Range 
to obtain aerodynamic data on the various fiseheties under consideration. 
These data were used in the preliminary development of the dispersion 


theory and are compared with the wind tunnel results in Appendix B. 
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Dynamic Supersonic Wind Tunnel Tests of 
Four Flechette Configurations 


i In order to obtain both static and dynamic wind tunnel data on 
flechette configurations, special tests were carried out at the 
University of Notre Dame which utilizes its unique vertical down flow 
supersonic wind tunnel and utilizes its one-degree-of-freedom pitching 
dynamic support instrumentation, Four flechette configurations were 
constructed and tested. The data from these dynamic tests was 
measured on a photo-comparator and reduced and fitted by using the 
Wobble program, The Notre Dame data on CM, and a CMs is 
in good agreement with the data obtained by the Ballistic Research 
Laboratory at small angles of attack and small mach numbers. At 
' the larger angles of attack, the Notre Dame data is as much as four 
times larger as the BRL data in damping und as much as two times 
larger than the Cy, data. Thus, the nonlinearities which have been 
uncovered in the dynamic wind tunnel tests are of considerable 
importance in evaluating flechette flight performance and in evaluating 


flechette accuracy and dispersion. No wind tunnel data was obvained on 


the important Magnus moment. This omission is considered extremely 


serious and it is recommended that future studies be carried out in 
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this area. It is also recornmended that the aerodynamic characteristics 
of the ditferent flechette designs be evaluated with a view towards 
improvement in performance and accuracy. 

Preliminary tests were carried out in obtaining the rolling 
motion of flechettes at the various angles of attack and in obtaining 
three-degrees-of-freedom wind tunnel tests where models were able 
to freely pitch, yaw, and roll. The exploratory rolling tests were 
carried out in the supersonic wind tunnel at Picatinny Arsenal, Good 
success was obtained on the basic configuration at small angles of 
attack. At the large angles of attack the sting support mechanism bent 
and thus had to be redesigned. These rolling tests have demonstrated 
that it will be quite possible to obtain excellent free rolling motion 
performance of flechettes at srnall and large anbles of attack using 
instrumentation at Picatinny Arsenal, 

‘Three ~degree-of-freedom dynamic wind tunnel tests were 
explored ina preliminary way in the Notre Dame vertical down super- 
sonic wind tunnel. In these tests the model was able to freely pitch 
and yaw and the afterbody with fins was able to roll freely. The fore- 


body however did not roll The tests were of marginal success but 


suggested that complete success could be achieved with more effort. 


It is specifically suggested that the new 3-D testing procedures originally 


explored at Notre Dame be continued in the Picatinny Arsenal and/or 


the BRL wind tunnels. 


It should be emphasized that the nonlinear aeroballistic dynamic 
stability coefficients obtained in the Notre Dame program represent a 
major finding which was extensively utilized in the performance analysis 
and accuracy computations. It is considered essential that all future 


flechette designs undergo complete dynamic wind tunnel testing and 
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range firings in order to permit accurate computations of the true 
dynamic flight performance, accuracy and dispersion of flechettes. 
Dispersion Theory of High Fineness Ratio, 
Cruciform Fin Bodies 
A complete jump and dispersion theory is setforth for the free 
flight performance of flechettes. The six-degree-of-freedom equations 
of motion are coded for various computer computations which indicated 
that the flechette accuracy theory accurately predicts the jump and 
4 
|, dispersion of flechettes. 
In order to determine realistic values for the initial conditions of 
flight and for the actual dispersion of flechettes, test firings are carried 
out in order to obtain special experimental data, The raw experimental 
" data is fitted by the least squares method and thereby placed into the 


form of initial flight conditions. These initial conditions are then 


applied to the theory. Six-degree-of-freedom numerical computations 


ment between the theory and test firing results indicate that the methods 


of data analysis and the flechette accuracy theory together provide a 


precise means of predicting the dispersion of flechettes. 


t 
| 

hy 

| are used to evaluate the dispersion of eight test rounds. The good agree- 
| 
| 

| 

| 

| 

| 

| 
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The analysis of the firing data indicates that the large initial 
conditions of flechette flight result from a strong impulse imparted to 
the flechette in the muzzle blast regime. It is found that if the transverse 
impulse imparted to the flechette is equal to an opposite angular impulse 
then the dispersion will be zero. Since these two impulses rarely balance 
and always exist, flechette dispersion is generally large. However, by 
controlling sabot design and muzzle blast, the transverse and angular 
momentums can be reduced and partially balanced thereby yielding 
excellent accuracy and low dispersion. 

Of particular importance is the invalidation of the classical 
maximum yaw theory long used in exterior ballistics. 

More specifically the complete jump and dispersion theory for 
flechettes has been reduced to three governing equations which represent 
flechettes having high, low and very low roll rates. These three theories 
were found to be accurate by evaluation against six~degrce-of-freedom, 
numerical computations of the equations of motion. It was found there- 
fore that they accurately predict the jump and dupsnsion of flechettes. 

The computer program undertaken to evaluate the flechette accuracy 
theory includes 201 special case runs carried out in four parts. The first 
part validates the theory with respect to the aerodynamic restoring and 
damping moments. The effect of these moments on dispersion was found 

to depend on the initial conditions. 


The second part validated the theory with respect to the aerodynamic 


Magnus force and moment. The effects on dispersion were found to be 


very small and of no consequence unless the total dispersion of a 


particular round was of the same order of magnitude as a Magnus 
effect. 

The third part validated the theory with respect to aerodynamic 
asymmetries (mass asymmetry, inertia asymmetry, etc.) and roll 
rate. All three theories were found to be validated in this phase and 
found to be quite accurate. Aerodynamic asymmetries causing a trim 
of 1° have little effect on the dispersion of fast rolling flechettes. 
Slower rolling flechettes were found to have in general increasingly 
large dispersion values as the roll rate decreased. It can be concluded 
that for flights which are prone to aerodynamic asymmetry and fin 
damage, a high roll rate is essential to low dispersion and increased 
accuracy. 

The fourth part validates the theory with respect to gravity. The 
theory indicates a lateral contribution to dispersion from gravity in 
addition to the vertical contribution. However, for the flechette this 
lateral contribution was found to be minimum. 

In general, the agreement between the flechette accuracy theory 
and the computer computations were cxccllcnt and account for the effect 
of the initial launching conditions as well as the static and dynamic 
stability coefficients and asymmetries. Further, simple equations are 


given in order to achieve the desired accuracy and optimization. 
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SUMMARY 


Ry an exploratory firing program, by a supersonic dynamic wind 
tunnel testing program, by the development of an accuracy theory for 
jump and dispersion, by computer computations and analys's, and by 
precision range firings at Frankford Arsenal, flechette accuracy and 
dispersion is explained, evaluated and improved. 

The firing program revealed the importance of fin and body damage, 
the blast region and saboting. The dynamic wind tunnel program yieided 
values for the important static and dynamic stability coefficients. The 
flechette accuracy theory was confirmed by numerical integration of the 
6-D equations on the high speed computer where the effects of initial 
conditions, stability coefficients and asymmetries was revealed and 
evaluated. Finally, by a flechette firing program in the new Frankford 
Arsenal Ballistics Range, excelent correlation between theory and 


experiment for flechette accuracy was obtained. 
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APPENDIX A 


EXPLORATORY FLECHETTE FIRING PROGRAM 


Two flechette firing programs were carried out at the University 
of Notre Dame. The first program was carried out in the Army Firing 
Range located under the football stanus in the Rockne Stadium, In these 
first firing tests the actual flechette and its sabot were fired at full 
hypersonic velocity using a mann barre] with sabot stripper. The firings 
were carried out with the assistance of technical personnel from Frankford 
Arsenal and under the direct supervision of Army ROTC personnel 
stationed on the campus and responsible for the Firing Range. These 
firings revealed two very important discoveries, By firing through light 
drawing paper yaw cards and by examining the impression left by the 
passage of the flechette, it was possible to obtain a positive confirmation 
that the fins were being seriously damaged and/or bent by the stripper. 
This finding was transmitted to the cognizant Frankford Arsenal personnel 
where suitable corrective changes were initiated and finalized thereby 
eliminat.ng the problem of fin damage. 

The second major finding of the first flechette firing program, insofar 
as university investigators were concerned, was the recognition of the 
tremendously intensive and long muzzle blast regime. While a standard 22 


projectile is fired in the range with little noise and little blast, the flechette 


system of basically the same weight but fired at large velocity yields a 
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tremendous concussion and a tongue of fire, blast and flame stretching 


some 3-4 feet. . The importance of the recognition of the strong blast 
region lies in its effect in disturbing the flechette at launch and thereby 
contributing to inaccuracy. 
The first firing program therefore revealed fin damage due to 
sabot strippers and a large blast region which contributed to jump and 
inaccuracy. 
The second flechette firing program carried out at the university 
| utilized an air gun in simple subsonic launchings. The range setup is 
| shown in Figure 1. The purpose of this special firing setup was to explore 
various saboting techniques. In this program sabots of both pusher design 
and puller designs were investigated. Also body inset sabot designs were studied, 
see Figure 2, Representative target data is illustrated in Figure 3 where 


1 effects of sabot designs are clearly evident. Various flechette and sabot 


designs are shown in Figures 4 and 5. 
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Figure 3. 
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Comparison of Pusher and Center Sabot Results Flechette Testing 
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Figure 4.) Producibility Sabot ~ R & D Plechette 


Various Puller and Pusher Sabets and 


Flechette Configurations 
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FOUR-FLECHETTE CONFIGURATIONS 


DYNAMIC SUPERSONIC WIND TUNNEL TESTING* 


| 
7 
| APPENDIX B 
a DYNAMIC SUPERSONIC WIND TUNNEL TESTS OF 
ABSTRACT 


| ; The linear values of the static pitching moment stability coefficient, 

| 7 CM and the damping moment stability coefficient, Cpy +CM 3) 

| are determined versus angle of attack for four flechette designs. The 
program is carried out in a vertical supersonic wind tunnel using a 

| one -degree-of-freedom dynamic testing technique. This method allows 
the model to go through free one-degree-of-freedom angular oscillations. 

Stability parameters are extracted from a film record of this motion 

f and the stability coefficients are computed using the WOBBLE computer 


program. Good repeatabilitv cf the results is shown for low angle of 


attack, 


*Prepared by Michael Garsik. 
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INTRODUCTION 


With the advent of more advanced analysis mmchnigues: today's 
aerodynamicist has the power to achieve a better understanding of the 
free flight performance of a flight vehicle. Data such as angular motion, 
jump angle and dispersion can now be extracted from free flight data and 
studied” so that previously undetected instabilities and design failures 
can be corrected. Obviously from this there arises a clear need for 


3,4 The random method used in 


development of free flight simulations. 

trying to solve the problems of stability and flight performance wouid 

prove dangerous and costly if full scale flight tests were conducted. It 

would be much cheaper and safer to experiment with new designs on 

models of the actual configuration. This presents the problems of simu- 
lating free flight motions so that data can be extracted and the new designs 
evaluated just as if the test were conducted cn a full scale model in free 
flight. 

Ballistic range firings was one of the initial attempts at a flight 
simulation technique. It involved taking photographs at various stations 
along a firing range of a model that had been launched from a gun. Be- 

cause of the limitations on the types of motion that could be observed, 
the lack of control of initial conditions, and other limiting factors, it 
soon became apparent that a more sephisticated method of simulation was 


necessary. Attention was turned to the wind tunnel. 


Attempts to study the angular motions of flight vehicles in the wind 
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tunnel began by mechanically reproducing them. This technique ran into 
several problems, in particular separating the driving mechanism 
response from the aerodynamic response and the fact that the technique 
is limited in that a mechanical response, rather then a free one, to the 
flow field is used. In recent years the most successful wind tunnel simu- 
lation technique, dynamic wind tunnel testing, has been developed. 
Actually there are several types of dynamic wind tunnel testing, The free 
flight angular oscillation inethod exhibits complete six-degree-of-freedom 
motion and needs no external support system, however certain limitations 
to this technique do exists. The duration of the simulation is restrictive 
hence length of the "flight" is very short. Also, a lack of control of initial 
conditions prohibits the study of particular flight modes. Another method, 
that of constrained angular oscillations, eliminates these disadvantages 
at the expense of introducing new ones. The most predominate disadvant- 
age is the interference effects of the support system on the response of 
the model to the flow field. This assumes that the problem of building an 
adequate support system can be solved. It is important to have control 
over the initial conditions and the length of the simulation run in order to 
simulate the free flight angular motions in the wind tunnel. Of course, the 
choice of which method to use depends on the careful consideration of the 
problem at hand and the experimentai limitations which could be allowed 
and not interfere with the test being carried out. 

With regard to the constrained angular oscillation technique and 


its use in the superso ‘c wind tunnel, several mcthods have been 


26 


developed. The gas bearing system is onc that is ideally suited to the 
study of low fineness ratio, non-finned bodies such as projectiles. It 
does not lend itself to the study of high fineness ratio finned bodies quite 
as well. One of the drawbacks of this technique is the high cost of 
construction and maintenancc of the system. The jewel bearing support 
system has been utilized in supersonic wind tunnel testing to observe the 
rolling motion of various models of flight vehicles. Such a system has 
been successfully employed in determining the roll damping moment and 
induced roll moment stability coefficients for different flight contigura- 
tions. 

This investigation is intended to determine the linear pitching mo- 
ment and damping moment stability coefficients of four flechette configu- 
rations in a supersonic regime. The study was conducted under a 
contract awarded to the Department of Acrospace and Mcchanical 
Engineering at the University of Notre Dame by Frankford Arsenal, 
Philadelphia, Pa. The contract deals with a .:tudy of the jump angle and 
dispersion of the flechette configurations, An underlying intent will be to 
document the constrained angular osciliation technique used in the 
supersonic wind tunnel tests. 

In order to study the performance of the flechette configurations 
and to be able to predict their flight path, a basic understanding of the 
stability of the rounds must be obtaine:I, Adequate stability prediction 
requires that techniques of flight simulation be used which will produce 


continuous results for supersonic conditions. 
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The actual steps taken in developing such a program of dynamic 


j 
wind tunnel tests were: 1) adapting the one-degree-of-freedom free | 
oscillation technique to the supersonic wind tunnel; 2) recording the one 
degree of freedom angular oscillations of the models in the supersonic 
wind tunnel by high speed photography techniques; 3) reducing the motion 
of the models to numerical values of angle of attack; 4) fitting the Aero- 
ballistic Theory to the angular data obtained to determine the stability 
parameters Kn, : ~NP? Pegi 5) computing the aerodynamic 
stability coefficients from Linear Theory using the stability parameters, 
model parameters, wind tunnel Mach number and density; 6) analyzing 
the interference of the support system by checking the repeatability of 
results, 

To accomplish the goals set down a unique method of supporting 
the pure pitch flechette models are utilized.’ It involved suspending the 
model in the test section of the University of Notre Dame's vertical 
supersonic wind tunnel and allowing it to go through free one-degree-of- 
freedom oscillations. The low friction in the system allowed continuous 
motions to be obtained and recorded and the stability coefficients to be 


extracted from the angular data. 


JN 
ix 


AEROBALLISTIC THEORY 


Axis Systems 

Two basic axis systems are used. The space fixed axis system 
(Figure 1) is the system in which the data is recorded. ‘The aeroballistic 
axis system (Figure 2) is the system in which the equations of motion 
are expressed. By choosing the x-axis of the space fixed system to 
coincide with the velocity vector the data is made directly compatible 
to the equations of motion, Frorn Figures 1 and 2 it is seen that @=@ 
and q= 6. Care must be taken in extending this comparison beyond this 
point. 

é 
The linear theory for a missile constrained at its center of gravity 


for one-degree-of-freedorm pure pitching is as follows. 


Linear Theory 


In the development of the Linear Thcory several assumptions are 


ie aS ae na eae 


made: 
1, Aerodynamic coefficients are constant dl 

sa 

2. Velocity and density are constant ! 

3. All angular motions except roll are small enough that the 


sinx = tanx = xX 


lm Lie ct io aor A a nei 


{ 

) 

f 

small angle approximations may be used: 
| 

cosx = 1 


4, The missile has mirror symmetry and trigonal or greater 


? rotational symmetry. 
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The fundamental differential equation of motion for the rotational 
motion is 
M = 16 (1) 
The sum of the acting aerodynamic moments, shown in Figure 3, which 
are assumed to vary linearly with angle of attack is 


M = M,@ +M a+ M,@ + Ms §. (2) 


q 


where 


= 
" 
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(3) 


Because of the selection of the particular axis systems and their 


orientation, Equation 1 can be rewritten as \ 
M = 1a (4) 


Equation 2 can be rewritien as 


M=M,a+ MytMe @ + Mg 8, (5) 
i 

7 

‘ Combining Equations 4 and 5 and rearranging 

* Ps M + M ee M 

: a - |_—a__& a &~ {S| a = My 8 (6) 


and 
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where 


M 
2 sale Ea | 
M; 5 


(7) 


(8) 


Solving for the homogeneous solution to Equation 7 assume a 


solution of the form 
a= Ke pt 


Differentiation of this yields 


ad = oKe Pt 


a = 7 Ke? 


(9) 


(10) 


Substitute Equation 9 and 10 into the homogeneous form of 


Equation 7 
s KePt 4 Ngee +N, ePF=0 
o7 +N, +N, =0 
1 2 


which has a solution of the form 


N 
ee ane | 1 Z. 
Sg eg pay PANG 


(11) 


For missiles in air the assumption that the products of stability 


2 
derivatives are negligible when compared to themselves (i.e. N, <<No) 


wre mane ot ne eee ate ce ee 


can be made. This is generally a good assumption and will be made 


here. Hence Equation 11 can be written 


Ny 


nee hyNy 


- + 


P12 


(12) 
PR EAE 
The homogeneous solution has the form 
(A ;+iw,)t Ay ~ iwg)t 
a= Ke + Ke (13) 


where 
CoM + Cry 
_ d 7) 1 

ae Mote | g @ Usd 21 se, 

] 2 

C usd} } 

w jase i (15) 
. 1 L 
scerras: (16) 
| 
| 


and 


‘ Ca, .) . 8B 1 
3 Ng * Ma) Suga ats) 


: Solving for the particular part of the solution of Equation 7 con- 


sider the steady state case of no pitching, Equation 7 would be 


N é (19) 
a=. Ma = Kg 


Ny M 


This is the particular part of the solution of Equation 7. The complete 


solution is 


(A; +iwy,)t (Ay + iwa)t 
a= Ke + Kye + K3 (20) 
where Ki and Ky are found from initial conditions and are 
we A / a - a K 
K. = %o 7 92,1% . $2,183 (21) 


1,2 : Fey meas 
91-2 ~ 21 (by 97% 1) 


Since the magnitudes of o1 and y are always equal and Kg is Constant 


l i K, = K 
Since <a 


| iwt: ae 
: e = cos wt + 1sinwt 


| Equation 19 can be written as 


a Ree 


whee. 


At 
: a= 2Ke cost + Ke 


a 


St. at | 


| : or in a more general form 


At 
a= Ke cos (wt+5) + K, (22) 


where 
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Equation 22 is the basic modal which will be used to fit the one- 
degree-of-freedom data. A physical representation of what Equation 22 
means and how it reduces the Tricyclic Theory to a pure pitching motion 
case is shown in Figure 4. The two arms K, and Ky have been replaced 
by’a single arm of length K where K=K +Ko. This arm is rotating at a 
rate w=w, and has an initial orientation of 6 =6 yj: Phe cosine function 
projects this arm onto the vertical axis of the aeroballistic axis system 
to give values of 8. This "projection" follows the pure pitch*ag of the 


model as if it would look when observed from the rear, 


Computation of Aerodynamic Stability Coefficients 
To fit Equation 22 to the angular oscillation data the WOBBLE 
computer program was used. This program fits the theory to short 
segments of the data in overlapping picces so that the stability para- 


meters Ny w) and Ky are determined as functions of time. 


Computation of Lincar Coefficients 


Using the velocity and model parameters (Appendix A) along with 
Ay: @®. and kK, the pitching moment stability cocfficient, CMa » and the 
damping moment stability coefficient, (C4 + CMs ), were computed. 


Equations 17 and 18 were used to compute these coefficients as functions 


of time. 
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Experimental Technique 


Four different configurations were tested, the Ground Point, Olin, 
Swaged Point, and Tracer. Schematic representations of the made] 


configurations are given in Figures 5, 6, 7, and 8 respectively. 


One -Degree-of-Freedom Wind Tunnel Test Procedures 

All of the tests were carried out in the University of Notre Dame's 
vertical supersonic wind tunnel shown in Figure 9. This wind tunnel 
features a vertical test section fitted with interchangeable stcel and 
glass walls. 4 steel wall was used on one side to give maximum support 
to the model support system and a glass wall was used on the other side 
to allow observation of the models. The basic idea behind the support 
system is shown in Figures 10, 11, 12 and 12a. 

To mount the model in the tunnel the following system was used. 
A length of piano wire 0.030" in diameter wag inserted through the hole 
in the glass wall and into a syringe tube. The purpose of the two syringe 
tubes, one on each side of the model, was to insure that the model would 
remain in the center of the wind tunnel test section after it was released 
and allowed to oscillate. After running the wire through the first syringe 
tube, it was pushed through a small hole 0.040" in diameter drilled per- 
pendicular to the longitudinal axis at the center of gravity of the model. 
The wire was then pushed through the second syringe tube and guided out 
of the test section through a hole in the steel wall. ‘he wire was secured 
outside the wind tunnel test section by a system shown in Figure il. On 
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figure 9. Vertical Supersonic Wind Tunnel 
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the glass wall side the wire was put through a hole in a flat washer 
placed flush against the giass wall. The hole in the washer had a dia- 
meter closer in size to the diameter of the wire then the hole in the glass 
wall. This cut down on the disturbance of the flow caused by the pre- 


sence of the hole in the wall of the test section. Next, the wire was 


secured by placing it through an Allen screw tube and a solder washer. 
An Allen screw tube is a long cylinder with a smal] diameter hole drilied 
along its longitudinal axis and three small holes drilled perpendicular to 
the axis. ‘These holes have been tapped to accommodate Allen screws 
Which can be tightened to clamp down on the piano wire and hold it in 
place. A solder washer is a short cylinder containing two small dia- 
meter holes, one at the center and one near the outer edge. After running 
the wire through the center hole it cin be bent around into the second 
hole at the outer edge. This hole also has a small hole drilled perpen- 
dicular to it and tapped to hold an Allen screw. As in the Allen screw 
tube, this Allen screw can be tightened down on the support wire to 
secure it. This system holds the support wire on the glass wall side of 
the test section. 

After running the wire through the steel wall side it was pushed 
through a flat washer identical to the one on the glass wall side. A 
tightening tube was placed next in position and the wire was guided 
through it. A tightcning tube is two concentric cylinders which are 


matched by threading. The length of the tube can be adjusted to the 


desired size by rotating the outer tube about the inner one. Finally, the 


piano wire was Suctiea using an Allen screw tube and a soldier washer. 
This completed the setting up of the support system. One advantage not 
already mentioned comes to light at this point. The models could be 
easily removed and inserted into the test section of the wind tunnel. 
, After the support system wés in place the tension in the wire was 
adjusted by changing the length of the tightening tube. The tension was 
set so that the model would not change its vertical position after the 
tunne! was turned on. The model was then rotated 180° and held in place 
by a retaining mechanism shown in Figure 13. This system consisted of 
| an extendable retainer which was placed around the ‘nose of the model to 
secure it at its initial angle of attack. The retainer was connected to a 
release wire which could be manually operated from outside the tunnel. 
When the release wire was pulled the retainer would slip off the nose of 
the model and the model would be free to oscillate. 

To record the oscillations of the model a Wollensak Fastex high 
speed motion picture camera was used. The camcra was set-up as shown 
in Figure 14 on the glass wall side of the test section. Two floodlights, 

: one just above the camera position and one above the inlet of the wind 


tunnel, were used to provide maximum lighting of the mode! in the test 


for three seconds with a lens opening of f5. 6. 


| section. The camera was operated at a speed of 3000 frames per second 
| Upon completing all preparations the tunnel was started following 
| 
( 


the procedure in Appendix B. The retainer was pulled back and the 


H 
: subsequent angular motion of the model was recorded. 
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One - Degree -of- Freedom Data Reduction Procedure 

The one-degree-of-frecdom oscillations were converted to 
nuumcrical values of angle of attack in the following manner. Two 
reference dots at a known distance apart had been placed on the steel 
wa}l in the test section behind the model. These dots were included in 
each frame of the film record of the angular motions of the model. The 
dots were placed such that a horizontal linc running between them was 
above the highest point that the model with the largest radius would 
reach. For each configuration the radius of oscillation, the distance 
from the pivot point to the nose of the model, was also known, The re- 
lative coordinates of the reference dots and the nose of the model were 


determined from the data film using an optical comparator shown in 


Figure 15. A computer program called REDUCE, presented in Appendix 


C, using these coordinates and the known conversion distance between 
the reference dots was then employcd to produce a time history of the 


angular oscillations of the model. A schematic of the reduction coordi- 


nates is presented in Figure 16. 


Velocity Determination Technique 


Since all the tests were not conducted on the same day it was 


necessary to determine the velocity in the wind tunnel on the particular 


day the test was conducted. A method of measuring the static pressure 
in the test section was necessary to do this. A system like the one in 


Figure 17 was used to do this 
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A long thin tube was placed in the inlet of the wind tunnel] and 
lowered until the end reached the support strut just below the test 
section. This end of the tube was secured to the strut to help maintain 
the position of the tube near the center of the wind tunnel test section. 

A small hole had been drilled in the side of the tube to coincide with 

the position of the model when the tube was in place. The end of the tube 
in the tunnel was sealed and the open end was connected to a manometer 
by a length of plastic tubing. This upper end was fastened so as to put 
tension on the tube and prevent it from moving about in the test section 
when the pressure measurement was being taken. Any movement of the 
tube would affect the pressure reading and produce an incorrect value 
of the velocity. 

Before starting the tunnel a tare reading was made on the mano- 
meter and the stagnation or total pressure was taken from a barometer. 
Since the manometer scale did not coincide ‘viih that of the barometer 
the tare reading and barometer reading, which should have been equal 
had their scaics coincided, were different, This dieconké was a 
correction factor which would have to be added to the pressure reading 
taken when the tunne} was on to give the actual static pressure, ‘The 

tunne! was turned on and the pressure was recorded. IJaviny ail] of these 


pressure readings the ratio of static to total pressure could be solved 


for using the following formula: 


Pp _ Pread + (Pt > Prare) 
Bos, i, ee 
Py Ye 
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: | : where 
a ; 
. (PU =P ) = correction factor —.--.- sates , 
i C tare | 
YS Once this ratio was known the Mach aumber could be found in the } 
: ; Isentropic Flow Tables of Reference 8. For ihe REDUCE computer , 
| program it was necessary to put the velocity in units of feet per second : 
t from the Mach number. A sample calculation of this is shown in Appendix , 
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ONE-DEGREE-OF-FREEDOM TEST RESULTS 


pt nr nmenpne-sepree—of-Ercedom Data Reduction 


The WOBBLE computer program was used to fit the Aeroballistic 
Theory to the one-degree-of-freedom angular data obtained from the 


REDUCE program. The data was fitted in segments of 1.8 cycles and 


the stability parameters Ki Kans Ay and Ww, were determined by 
| WOBBLE at a time interval of 0.03 seconds. The stability parameter 
Ky , the trim arm, is analogous to the K3 in the Linear Theory and is 
due to aerodynamic asymmetries in the configurations. The average 
pexcent error of the fitting of the thcory to the data for all the tests 
carried was less then 3%. A representative plot of probable error (P.E.) 
versus time is given in Figure 18. The stability parameters were ob- 
tained from the fits as functions of time. Plots ef the stability para- 
| meters versus time for all the configurations are presented in Figures 


19 through 34. 


ae One -Degree-of-Freedom Stability Coefficients 


The pitching moment and damping moment stability coefficicuts, 


~~ 
a ed 


CM and (Cm, + CM my »were obtained versus tirue from the WOBBLE 
fits. Plots of the mean valucs of the coefficients per fit versus mean 
angle of attack per fit are presented in Figures 35 through 42 for all the 
configurations. These plots given an approximation of how the coeffi- 
cicnts vary with angle of attack. Included on these graphs are Ballistic 


n 
t 

i 

? 

Ey oA Range La! oratory (BRL) results for the respective configurations and 
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coefficients. The BRL. data was plotted at an angle of attack of 2°, 

Figures 43 through 50 are plots of BRL results for the stability coeffi- 
cients versus Mach number for all the configurations, Mean values of 
the Notre Dame results for the coefficients at low angies of attack are 
included on these graphs. The Notre Dame data was plotted at a Mach 


number of 1,3 which was an average value of all the tests carried out, 
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An important point which should be brought up at this point is the 
discrepancy in the definition of the damping moment stability coefficient 


between the two sects of results. The computations in this investigation 


| 


were carried out using a factor of ( 5. -) in the definition of the damping 
moment stability coefficient (see ee 3). The BRL definition used 
a factor of Gey) Causing the respective computed values of (CMg*CM a) 


to be off by a factor of two. To account for this and allow the results to 


be directly compared, the BRI, values of (CM gt CM y) were increased 
by a tactor 2 before plotting, This essentially gave all the velues pre- 


sented a uniform definition and allowed the comparisons of the results 


to be made, 
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CONCLUSIONS 


Single-degree -vi-freedom dynamic supersonic wind tunnel tests of 
four flechette configurations has been presented. Linear values of CMe 


and (Cyy_+ Cry.) were determined from stability parameters acquired 
Mg Mg P 


from the data taken during the tests. These values of CMa and 
CMgt CM ¢, showed good repeatability and were compared to results from 
the Ballistics Range Laboratory (BRL) for the same designs at low angles 


of attack. Over the range of comparison the agreement between the two sets 


of data was shown to be quite good. 

The repeatability of the results was a good indication of the absence 
of frictional effects and interference effects to the flow which might have 
been caused by the support system. In Reference 7 it was shown that this 
excellent one-degree-of-freedom dynamic testing technique can be easily 
extended to include the determination of nonlinear values of the static 
pitching and damping stability coefficients, Also, because the model is 
suspended at its center of gravity, there is no reason why this technique I 


could not be moved into a horizontal supersonic wiud tunnel if necessary. 


APPENDIX A 


MODEL PARAMETERS 


Diameter = .012 ft. 
t Le . 000004647 slugs-ft2 
Mass = .0003647 slugs 


| 
3 Grount Point 
Radius = 1.69 in. 


i Olin 


Diameter = .O119 ft. 
Mass = .000303 slugs 
I, = .000007040 slug-tt2 
Radius = doz ing 


Tracer 
. Diameter = .01533 ft. 


Mass = .0004520 slugs 9 
f I, = .0000105570 slugs-ft 
: Radius = 2.11 in. 


; Swaged Point 


. 
; Diameter = .012 ft 
Mass = .0003933 slugs 
\ ly = ,000006041 siugs-ft 
f Radius = 1,69 in. 
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APPENDIX B 


SUPERSONIC WIND TUNNEL OPERATING PROCEDURE 


Starting Procedure 

1. Open valve to compressor manifold for wind tunnel to be used - 
make sure that the other wind tunnels are either shut off or blocked from 
the manifold. 

2. Inform University Power Plant of intention to run compressors. 

3, Turn cooling water on (one valve near wall inside laboratory). 
4. Turn each compressor shaft to make sure they are free to 
rotate, 
5, Check oil level for cach compressor (oil level should be above 
gear). 

6. Check oil pump for each compressor i.e. depress six plungers 
and observe oil bubbics. 

7. Adda few squirts of No. 5! oil to nole in top of shaft bushing. 

8. Check mercury manometer tubing in compressor room to make 
sure it is connected, 

9. Turn master power switch on for cach compressor. 

10, Start one compressor - allow at least one minute after com- 
pressor comes up to speed before starting the second compressor and 
allow another one minute after this compressor comes up to speed before 
starting third compressor, 

lil.. If mercury manometer reads more than 18 inches, Shut Down 


Immediately, 
96 
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Shut Down Procedure 


1, Shut compressors off one at a time at one minute intervals, 
2. Turn master power switch off for each compressor. 


3. Shut compressor cooling water off. 


4. Inform University Power Plant that compressors have been 


shut off. 
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APPENDK C 
ONE -DEGREE-OF-FREEDOM ‘CEST RESULTS 
ON R&D FLECHETTES 
The model* was initially disturbed to an angle of attack of approxi- 
mately 180° and then allowed to oscillate freely. The resulting angular 


r 


motions were then recorded by a high speed camera technique. 


One-Degree -of-l'reedom Data Reduction 

The "Wobble" computer program was used to fit the one-degree-of- 
freedom Acroballistic Theory to the angular oscillations obtained from the 
moving camera technique. This data was fitted in segments of 2.2 cycles 
with each segment containing approximately 25 points. The stability para- 
meters K,, Kins Ay W,, were determined by the Wobble program at a time 
interval of 0.015 seconds. 'The average percent error of the theory to the 
data showed anerror of less than 3%. A representative plot of probable 
error of fit vs time is shown in Fig, ia. 

The stability parameters were obtained from the fits as functions of 
time, representative angular oscillations, probable errors of fit, and 
stability parameters are presented in Figs. 2 through 6. The resulting 


stability coefficients versus time are presented in Figs. 7 and 8. 


One~Degree -of-Freedom Nonlinear Stability Coefficients 


To get an indication of the nonlinearity of the stability coefficients . 


— *Figure lS 
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with angle of attack, the one-degree-of-freedom Nontineay Aeroballistic 

ng Theory was employed. Using this nonlinear theory, the stability coefficients 
were determined as polynomial functions of the angle of attack. Repre- 
sentative plots of runs made are presented in Figs. 9 and 10. 

Both C,, 4» the pitching moment coefficient and Cmng tCig the 

damping moment coefficient were found to vary nonlinearly with angle of 

attack, Beth were found to be highly repeatable. om = varied no more 

than 2%, about its mean while Cig + Cand varied less than 5% about 


its mean. 
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APPENDIX C 


DISPERSION THEORY OF HIGH FINENESS RATIO, 
CRUCIFORM FIN BODIES * 


A complete Jump and Dispersion Theory is developed for free 
flight vehicles. Six-degree-of-freedom computer computations 
indicates that the theory accurately predicts the jump and dispersion 
of flechettes. 

The initial conditions and dispersion values are established by 
range test firings, The raw data is fitted by least squares method and 
put into initial condition form. Initial conditions are applied to the 
theory and 6-D numerical computations to evaluate dispersion for eight 
test rounds. The results are compared to test firing target data. The 
agreement between the theory and test results indicate the data analysis 
and theory provide an accurate means of predicting dispersion of flechettes. 
Analysis of the firing data indicates that the initial conditions result from 
an impulse imparted to the flechette in the muzzle blast. The transverse 
impulse imparted to the flechette initially must be equal to the angular 
impulse to obtain zero dispersion. Other disturbances in the blast region 
such as sabot separation influence the initial conditions and hence 


dispersion, First inaximum yaw theory is discussed and disproved. 


*Prepared by I.awrence E. Lijewski. 
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INTRODUCTION 


The accuracy and dispersion of free flight vehicles has been a 
problem in aerodynamics and ballistics for many years. Until the present 
time, the primary investigations into causes and effects of jump (the angle 
between the line of boresight and the line connecting the point of launch 
with the instantaneous position on the trajectory, ) and dispersion have been 
directed toward projectiles and, in particular, artillery rounds. A full 
program to investigate jump and dispersion characteristics of low trajectory 
iinned bodies has been lacking and therefore is the subject of this disserta- 
tion. The purpose of this analysis is to develop a basic understanding of 
the parameters causing the jump and dispersion of flechettes. The flechette, 
being a gun launched finned body, requires a different approach to the 
problem. The old concept employed in the analysis of the dispersion of 
artillery rounds is that the dispersion results from initial launch disturbances 
imparted by the gun to the ahelle?* This concept is no longer valid for 
flechettes since the flechette is a fin missile, sabot launched, and its 
dispersion must be tied to the disturbances it encounters when clearing the 
muzzle blast and sabot separation region. In addition, asymmetries are 
more prevalent in finned bodies than projectiles and a finned body is more 
apt to be influenced by the blast. These factors must be taken into account 
by a theory involving finned bodies. 

In order to develop this new approach, (1) a theoretical expression 


for jump and dispersion had to be developed, (2} the theory had to be 
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validated, (3) free flight feet firings had to be undertaken and initial 
condition data extracted, and (4) the test firing results had to be 
correlated with the validated theory. The Jump and Dispersion Theory 
was developed, in general, for both fin and spin stabilized missiles in air. 
The theory -ncludes the effects of: initial conditions, Magnus, aerodynamic 
asvinmetries, and gravity. In the past, theory development for projectiles 
included only initial angle of attack and initial angular rate!+3 Initial 
trausverse velocity was considered non-existent@ or negligible. Zaroodny” 
included a linear momentum term to account for any transverse motion of 
the projectile but attributed it to the gun during recoil, Any transverse 
impulse imparted to the projectile by the blast was ignored. Other authors, 
including Sterné attributed the jump only to bore clearance and therefore 
: only included, effectively, the initial angle of attack. Magnus effects were 
always neglected in previous studies cither due to lack of familiarity with 
| the subject or lack of data. In general, all cross-forces, except lift, were 
neglected mainly for convenience sake, Zaroodny, however, cautioning 


against wholesale simplifving said “it would seern desirable that cur 
ving 


formulas allow us to include these other forces as the experimental 


information on these forces becomes available." Aerodynamic asymmetries 


were neglected for projectiles but included in Murphy's work.© It was not 


25859 


until Nicolaides that all four factors affecting dispersion; initial 


angle of attack, initial angular rate, initial transverse position and ve locity, were 


put into one theory. The work presented here expands the work of Nicolaides 


to include all parameters affecting dispersion in detail. Three separate if 
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! equations comprise the theory to include the complete range of roll rates. 
Before, only high roll rates were considered; with the study of finned 
bodies, the roll rate range extends down to zero roll and accurate 
theories had to be deduced from known aerodynamic equations. 

To validate the theory,a six-degree-of-freedom trajectory computer 
program numerically integrating the equations of motion was utilized’? ee 
The validation consisted of four phases. The prccedure began with the 
most basic theory equation and consecutively added terms to validate the 
entire theory. Initial conditions, magnus, asymmetries and gravity were 
successively validated with roll rate and velocity varied in cach phase, 
Before the advent of adequate photographic material, obtaining test 
data was often difficult, At first, jump target data was taken separate 
from yaw data. The thinking was that the yaw data was part of the pro- 
jectile’s characteristics and not affecting jump. As photographic methods 

| improved, and theories developed, the data was correlated. The correla~ 

tion of the data was often a problem. A fit of the motion to a least squares 

micthod was difficult. Fowler, Kent, and Ilitchcock developed 2 method 

i that would plot the magnitude of the yaw separately from the orientation 

aud then fit the curves separately. A better method was developed by 


| McShane -Charters-Turetsky approximated the yawing motion to a circle, 


For projectiles the method has been refined and is an excelent method. 
| However, for finned bodies with not, always circular angular motions, a 
different method of data analysis had to be devised. Utilizing the free , 1 
flight data taken by test engincers at Frankford Arsenal on a number of 


4 
i 
flechettes, the least squares method was employed to fit the data pre- } 
ve 
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sented here. The nearly planar oscillations of the fleciette in the first 


13,14 and the 


few feet downrange were fit to a pure pitching motion 
‘position downrange fit to a thira order polynomial. From these results, 
angle of attack, angular rate and transverse position and velocity were 
détermined for the first few feet downrange. Before, there was some 
controversy as to whether or not the least squares fit could be extra- 
polated back to the muzzle. Zaroodny contended that the x =0 position 

had to be taken out of the blast region to allow the aerodynamic equations 
to be valid. On the other hand, Kent, Hitchcock, Fowler and Sterne held 
to the fact that the free flight region began the instant the projectile left 
the bore. In the analysis of flechettes the position x=0 is taken somewhere 
downrange aicer the sabot separation sequence has occurred, This is seen 
to be 3to 5 feet downrange and assumed clear of any muzzle blast effects. 
: The striking shortcoming of previous works is the lack of correlation 
between test data and valid theory. For the flechette, correlation between 
the theory and test data was undertaken as well as correlation between 

test results and first maximum yaw data. Currently, the first maximum 

i ya‘ theory! is held by some to be an accurate method of predicting 
dispersion. This theory disallows any influence of initial angular rate, 
transverse position or velocity on dispersion. The dispersion analysis 
presented here disproves this theory with actual test data. The details 


of each of these aspects of this program are developed in the following 


sections, 


ea 


DISPERSION THEORY 


Dispersion relationships for free flight vehicles are embedded in 
the trajectory equation of any such aeroballistic body. To evaluate the 
trajectory equation and thus the dispersion, the linear second-order 


differentia] equation of angular motion is a logical starting point. 


= oa = o ipt. = 
w a NW Nowe NLS +N, (1) 


where Ni): No, Na, and N, are constants. 
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In this discussion of dispersion theory, it is assumed that, 


(1) total velocity, Ug» is constant, equal to u in the’ 
theory development. 

(2) all force and moment coefficients dependent on angle 
of attack are consiceréd to be linear with angle of attack. | 

(3) all force and moment coefficients independent of angle . 


of attack are considered to be constant. 


i 
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(4) a linear relationship exists between x (distance down 
range) and time for the non-drag case. 


(5) roll rate, p , is considered to be constant. 


(6) products of force and moment derivatives are negligible, 
except those involving Z, and M, . 
€ € 


=| 
| Utilizing these assumptions, and the binomial expansion of (Z,,,-m) K 


2, 3, 4 and § become: 


a er eee eset teal ei eel oil 


Z,, + ipZ M+ u ipl 
: woo ‘ ere ; i (2a) 
| y y 
+ipM j Z,, + ipZ 
Ny ® > u Ma RM py + be wt 'P “pv (3a) 
if I m 
kg y 
| »- ipZ,.6 UMg 5, 
Ni, = ele a (4a) , 
y y 
! 
\ i 
y 
7 *% t ‘ 
7 The solution to Equation 1 is that of tricyclic motion; that is, 
i Pe EE Fal t. = t,.trelltl_.F 
: We Ke + Ke gat + Ke + Ky (6) 


where the complex coefficients are: 
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The trajectory equation for free-flight motion: 


- 


and 


= (# - iuq) (11) 


An expression for q is obtained from the equations of motion 


a) - -™Mm <~ i i 5 e 
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a (123 
- Z+Z- Zu +ipZ., i7,. 6 . 
qe -iw E a) +iw pee =| |e | alpt  ¥ 
m mu mu | u 
yielding a solution of the form: 
. - . : + = % eG fy : 
Ss ke + koe 4 k,€ Kt + kt + Ke (13) 


where the entire 2xpression for ihe solution is: 
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The term G. is of an order of magnitude 10°? and thus is 


veglected from all further discussion. This reduces 14 to: 
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a Ry further inspection, terms with g? and ¢5 will be negligible since 
Q BOM 


they coutain products of force and moment derivatives. Equation 15 


becomes: 
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Uquation 16 contains only the significant terms in dispersion theory. 


aa This equation is valid for all values of roll rate 
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| High Roll Rate Theory 


For roll rates greater than 100 rad/sec, Equation 16 reduces to 


an approximate solution. Integration of the double integral gives: 


t ¢t 
? ipt 
elf drdr = fo . t -_t (17) 
Jf Gp? Pip)? 


For high roll rates, the first and third terms go to zero, leaving only 
the second term to affect dispersion. Applying this approximation to 


Equation 16 
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! where, by applying previous aerodynamic relationships: 
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Substituting 19 and 20 inte 18 and expanding the various terms: 
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Employing assumption 6 , 
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The mil-relation offers a method to define the Jump Angle from Equation 


22, 
Jump Angle = = (103) (23) 
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Equation 24 gives an approximation for the Jump Angle for high roll rate 
cases with gravity, at any position x down range. 
J.ow Roll Rate Theory 
Kor roll rates less than 100 rad/sec but having a parameter, pt, 
greater than 1, Equation 16 can be reduced to another approximation. 


As before, integration of the double integral yields Equation 17 
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For low roll rates all three terms are significant to dispersion. 


Equation 16 now becomes: 
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arm, or rolling trim vector must be separately examined, 
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From Equation 8, 
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or 
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(ip)* ~ ip($,+$2)+o1¢2 
Numerical inspection of the three denominator terms indicates that 
the first two terms can be neglected. Each term is not only less than 1% 
of the third term but also they're subtracted from one another to make 


their contribution even more minimal. Thus Ky is approximated by, 


3 md Sar he nde ee ea a 
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for low roll rates, the second term in the numerator and the first term 


in the denominator dominate all other terms and become the only signifi- 


cant terms, Thus, 


; (26) 


approximation holds true for applicable terms in 


! Equation 25 , thus reducing the jump angle equation to: 
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Combining terms and dropping the negligible second last term, 
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Equation 28 accurately approximates the jump angle for roll rates: 


p> 100 rad/sec 


pt > 1.0 


Very Slow Roil Rate Theory 
For very low roll rates; that is, p> O and pt<1', Equation 28 is 
again applicable. 


One approximation is used, however, and that is that cos CE) 
x ee: 
and sin( B) arc approximated by power series. 


2 
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Substituting‘and simplifying, 
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VALIDATION OF THEORY 


The theoretical expressions for Jump Angle; Equations 24, 28 and 
30; show that the dispersion depends on the initial conditions, aerodynamic 
coefficients, distance downrange , and mass parameters. Dispersion tor 
this theoretical analysis is defined to be the deviation from the line of fire. 
By analyzing only one Flechette configuration to validate the theory, the 
producibility Ground Point, and taking all cases to be evaluated at 1000 
feet downrange , then the expression for the Jump Angle can only be 
aftected by the initial conditions and acrodynamic coefficients. 

To assure that the three equations for Jump Angle are valid and to 
show the effects for various initial conditions and aerodynamic coefficients, 
the expressions for the Jump Angle were evaluated for a series of cases 
and compared to numerical integration of the six-degree -of-freedom 
equations of motion, (6-D). A sample case run can be found in Appendix 
A-2,. ‘The series of cases is broken down into various phases of develop- 
ment. Phase I considers various initial conditions but with only the 
restoring and damping aerodynamic cocfficients. This phase validates 
ihe use of initial conditions alone. Phase II utilizes a set of constant initial 

conditions, except for roll rate, and constant restoring and damping 
coefficients, while varying Magnus coefficients to determine their 
influence. Phase [I brings into consideration all the acrodynamic 
coefficients to include the configurational asymmetry coefficients. 


Different coefficients are used by varying the initial velocity and roll 


144 


manny a vee one me 


SS 


H 


rafes are varied to evaluate high, low, and very low roll theories. 
Phase IV considers the effects of gravity for various initial velocities 
and roll rates. No configurational asymmetries are used in order to 
isolate the gravitational influence. Values for ali coefficieits are 

found in Appendix Al, as well as other data including mass parameters. 
Since computations were done at 1000 ft downrange, the Jump Angle in 
mils is equivalent to the deviation from the line of fire in feet for all 
presented cases. The axis system used throughout this analysis is 


illustrated in the list of symbols. 


Phase I 


To validate the effects of initial conditions with restoring and 
damping coefficients only, 36 cases were evaluated using the high roll 
rate theory, liquation 24. The cases are divided into 4 sections 


isolating different initial conditions and their effects, 


Cases 1-9 


The first section shows the effects of roll rate and velocity with 
_ 


7 —> 
—_ e 
zero So» a, and QL 
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TABLE I 


THEORY VALIDATION, RESTORING AND 
DAMPING MOMENTS, CASES 1-9 


Coefficients eo 


nitia] Conditions 


: | 


ee 
5000 
0 


31416 
5 18850 Al 
6 
7 Poa RS 0 
1000 


8 QO | 18850 
Ee 0 6283 


Table I clearly indicates that no deviation from the line of fire occurs 


mn>aq 


if S.; Qo» and a, are set to zero. Roll rate and velocity changes have no 
effect on the Jump Angle for this particular situation. This is a trivial 


solution, it being obvious from inspection of Equation 24. 


Cases 10-18 


ye second section gives the effects of initial translational velocity, 
—_ 


Tl 
So = y + iz, with various roll rates and velocities. To assure the solution 
is correct in three dimensional space, the initial transiation velocity is 
given in both y and iz directions. Equation 24 reduces to: 


ae. o £000: = 
JN coo ie 
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THEORY VALIDATICN, RESTORING AND 


) 
TABLE If : 
DAMPING MOMENTS, CASES 10-18 


ote oe Coefficients 
Initial Conditions 


a 


5 100 is 
00-15 Fo | 1aa50|so0 
100 ‘ 
(fp [oo 
100 
‘W100. 2544 }100.000 +} 
100. 764 i]100.000 i 
| ) 100. 257 4{100.000 + 
100. 765 i}100.000i 
erates! LAESeneer ee | 
! The correlation between the theory and the 6-D integration for 
| Cases 10-18 is excellent as shown in Table II. The Jump Angle is scen to 
, Le affected by velocity but not roll rate, as would be expected from the 
| i : i 
| reduced Jump Angle equation. Figure | illustrates the deviation frorn the 


line of fire for initial velocities of 5000 ft/sec (Cases 10-12), 3000 ft/sec 


| (Cases 13-15) and 1000 ft/sec (Cases 16-18). Since the theory and 6-D 
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igure 1. Dispersion: Phase I Cases 10-18 


a 


i 148 ; 


100 50 


TDM PE ware nue mee 


120 


th) 800 10090 


1 - CASES 10,11, 12 
2. = CASES 13,145,195 
3 ~ CASES 16,17,18 


Ma2zr- AMODD- 


X~FEET 
200 400 600 800 1900 
SSS: 
ee ee 1 


f 
Maer AM 


Figuie 2. Trajectories, Cases 10-18 
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are so close, they are plotted as one point. Figure 2 illustrates the 
trajectory in both the x-y and x-z planes. The deviation from the line of 
fire. is linear with distance downrange in both planes. This would be 


expected with no gravitational force acting. 


Ree es wy me 


ree 


Cases 19-27 
The third section gives the effects of initial angle of acack, ae» with 
various roli:ates and velocities. Again a complex initial condition is 
used to validate the theory in three dimensiona] snace, [quation 24 


reduces to: 


JA. = iat phe pcs eee 1000 
©\mud C Ply eo 
My”! mao YL oy 


Table Ill shows the range of error between the 6-D computation and 
the theory to be 0.036 to 6.040 miis in the y-direction and 0.05% tc 0.041 
mils in the z-direction. Although the y~-dircction deviations differ in sign, 
the error between them is approximately 0.00225 degress, an extremely 
small angle. This angle will give an approximate deviation of 0.04 feet 

~— 

from the line of fire at 1000 feet downrange, With the J. A. being so close 
to zero it can be expected that the signs may diu ‘rt due to computetiona! 
errors. ‘The results do show Jump Angle variance with both roll rate and 


velocity. The largest changes occur.as vclocity goes to LOGO ft/sec. 
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TABLE lil 


THEORY VALIDATION, RESTORING AND 
DAMPING MOMENTS, CASES 19-27 


TA. (mils) 


Theory 


-0.027 
+0.027i | 


~0.017 
+0.017i 


-0.006 
0.008 | 
-0.026 
+0.016i 
0.0334 | -0.005 
ee 
0.0374 | -0.073 
O. Ta | 40,073: 
-(). OO84 | -0O.044 
0.0214) -O.015 
0.053i | +0.015i | 


0.023+ 
0,058: 


0.0124 
0. 067+ 


Cases 28-36 
a 


The fourth section gives the effects of initial angular rate 3 


with 
varying roll rate and velocity. An angular rate of 250 rad/sec is used in 
both directions of the complex plane to test validity in three dimensional 


space. Equation 24 reduces to: 
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TABLE {V 
THEORY VALIDATION, RESTORING AND 
DAMPING MOMENTS, CASES 28-36 


Cocfficients 
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Figure 3. 
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Dispersion: Phase I Cases 28-36 
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Table IV indicates excellent agreernent between the theory and 6-D 


computations. Roll rate is found not to affect the Jump Angle appreciably 
but velocity does, as would be expected from the reduced Jump Angle 
equation. Figure 3 shows the dispersion pattern while Figure 4 illustrates 
the trajectories. Cases 28,29, and 30 are plotted as one point due to the 
small difference between them. Cases 31, 32, 33 and 24, 35, and 36 are 


plotted similarly. 


Phase II 
To validate the effect of Magnus Forccs and Moments on the 

dispersion of flechettes, 21 Cases were run varying the initial roll rare 
and Magnus Coefficients. All other conditions were held constant. The 
variance of Magnus coefficients with Mach number had to be chosen 
since no data was available. Arbitrarily, the ratio of Cag /CMpp was 
chosen to be the same as that of C, o/—M y* The Magnus Coefficients 
uscd are presented as functions of Mach Number in Appendix Al with 


\ only the values at Mach 4.5 tabulated here for identification sake: 


TABLE V 


MAGNUS COEFFICIENTS, 
AT MACH 4.5 
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Equation 24 now becomes: 
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Initial conditions used in this scction are consistent with those of other 
sections to provide a basis for comparison. Three cases of zero Magnus 
were run, One at each roll rate to provide a standard to judge the influence 
of Magnus. 

The effects of Magnus coefficients on dispersion are minimal as 
seen in Table VI. The variance between the zero Magnus cases and any 
other case is found not to be greater then 0.209 mils (or feet at 1000 
feet of range). In order to obtain the maximum Magnus effects, the 
largest possible Magnus coefficients were used. Hence, Cy = 34.8 

PB 
and “Mpg = 110.0 arc the largest possible COCPREIC NES since cases 40 
and 49 become unstable. Table VI indicates the effects (for positive 
Magnus coefficients) 
(1) increasing horizontal dispersion with increasing p 
(2) decreasing vertical dispersion with increasing p 
(3) increasing horizontal dispersion with increasing Magnus 


(4) decreasing vertical dispersion with increasing Magnus 
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TABLE VI 


THEORY VALIDATION, MAGNUS, CASES 37-57 


a 


mMur>g 


9944 Pa 


.042i-Theory 


18.003+ 


18, 0821-Theory 
7 1T.3105 
17.944: | 
18. 141+ 
17.9031 
| 17.909+ 
17.942: 


18.1234 
17.913i 


18. 063i 


L7.954i 


17, 826+ 17. 910+ 17,9774 

18.2331 18. 144i 18. 0594 
17.8994 |, . 8204 
17.954i | 17, 942i 17. 931i 


157 


\y 


ne ee eee 


el a ron 


Z-FEtT 


18.2 


Figure 5. 


CmMéps8=4+90.0 


HEORETICAL,CASES 46,55 
HEORETICAL, CASES 47,585 
HEGRETICAL, CASES 48,57 
-D COMPUTATIONS 


1 
+ 
7 
6 


xec , 
ris 
P-INCREASING 


ere 


Dispersion: Phase II Cases 46,47,48,55,56,57 
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Figure 6. Dispersion: Phase LI Cases 38,41,44,47,50,53,56 


| (for negative Magnus coefficients) 


(5) decreasing horizontal dispersion with increasing p 

(6) increasing vertical dispersion with increasing p 

(7) decreasing horizontal dispersion with decreasing Magnus 

(8) increasing vertical dispersion with decreasing Magnus 
For example, Figure 5 illustrates the effects of roll rate for constant 

; Magnus coefficients of + 90° (1,2,5,6 above). Figure 6 illustrates the 

effects of Magnus for a constant sample roll rate (3,4,7,8 above). 
Obviously, when only a 0.209 mil maximum deviation due to Magnus 
occurs when the situation is geared toward finding the largest effect due 
to Magnus, smaller deviations due to Magnus would be found in actuai 
situations. It can be concluded that Magnus has no large effect on dis- 
persion although it could be significant if the total dispersion is close to 


. Zero. 


Phase ITT 
To validate the effects of aerodynamic asymmetries on dispersion 


of flechettes, a large number of cases were run varying roll rate, velocity, 


The asymmetries coefficients were selected to allow 1° of non- rolling 


| trim to exist while the flechette was in flight. Vhe asymmetry coefficients , 


Cyt Cor, Cup and Cyn are presented in Appendix A-1 as a function 


ler of Mach number. ‘The variance with Mach number was chosen arbitrarily: 


the ratio of asymmetry force to asymmerry moment identical to the ratio 
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of Cz, to CM 4 The wide range of roll rates makes mandatory use of 


all three dispersion theories. The governing equations are presented as 


they apply. 


Cases 58-90 
The first set of cases utilizes zero initial disturbances while vary- 
ing velocity and roll rate. For roll rates of 31416 rad/sec down to 100 


rad/sec the High Roll Rate Theory yields the governing equation, 


ie 
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For roll rates: p< 100 rad/sec and pt > 1.0, the Low Roll Rate Theory 


takes effect: 


Finally, the very Slow Roll Rate Theory applies for values of pt< 1.0: 
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Tables 7, 8, and 9 list Cases 58-90: 


161 


, 
Sidepenuiiateiehsstnasiathddecinadiecieiechata ernie ee er Oe tren Sin nancmeeacsianibaal LR oS teaad ath els nisbilbitinemantahe 
an ence i an en nt cn cea p/n ae ae 


Armaan tnt are hk ema ee Panne mae 
pa tans — Acie sates eemmmapmaeiaey .rinmenn 


. 
ee 


TABLE Vil 


THEORY VALIDATION, ASYMMETRIES, 
CASES 58-68 


Coefficients \ 
Initial Conditions J.A. (mils) 


Theory 


0.018- 
0.013: 


0.030- 
_0.0271 
0.060-| 0, 064- 
0.1271] 0, 130i 
0.997-| 1,013- 
0.9921} 1.009% 
1,620-| 1. 688= 
300 I.721i} 1.683% 


4.574-| 4.675- 
) 1000/5000 Al At | ALI a’ Boel 4.9754 
0 1 A 
| 


§.666-] 8. 780- 

ee 12, 489i 

|} 20, 669-| 21, 150- 
26.418i| 26.9271 
7.973 | -8.210 
-62.197i| -63. 210i 


-49, 706i | -50. 4271 | 
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TABLE Vili a 


THEORY VALIDATION, ASYMMETRIES, 
CASES 69-79 


<= ——— 
C Coefficients 
K Initial Conditions 2 
Cr 
S Ze 
E 1 


0.008 - 
0, 004i 


0. 009- 


6 0. 004i 
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70 
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1.9844 
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%,217 
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1.029" 1-42,27 33 
-41. 985i 


-33.014 
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TABLE IX 


THEORY VALIDATION, ASYMMETRIES, 
CASES 80-90 


TN (miis) 


Theory 


Unstable 
18850 Unstable 
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O.QLOi 
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Evident from Tables VII, VIII, IX is the fact that roll rate has 
tremendous influence on the dispersion of flechettes with aerodynamic 
asymmetries. Figures 7, 8, and 9 illustrate the dispersion pattern for 
these cases. The 6-D computations and theory are in very good agreement 
considering the large deviations invoived. It should be noted that the actual 
flechette with its velocity approaching 5000 ft/sec is affected very little 
by aerodynamic asymmetries. However, if the flechette were only t. roll 
very slowly, large dispersion ranges in excess of 60 mils could occur. 
Velocity also has a noticeable effect an dispersion. Figure 10 shows the 
three theory curves fim iigures 7, 8, 9 in composite to illustrate 
velocity effects. A saraple trajectory, Case 79, is shown in Figure 11, 
illustrating the curved path of flight. ‘This is typical of trajectories 


involving acrodynamic asyrnmetries. 


Cases 91-122 
Yo show the relation between the effects on dispersion for initial 
transverse velocity and acrodynamic asymmetries a second sct of cases 


were run, Roll rate and vejocity were varied as in the first set of cases, 
—. 


— 
2 3 7 — ° jae 2 
but 55 was set at (100 4 lOQOi) ft/sec with a, = 0 and o=0. Tables X, XT, 


and XII list the results. For high roll rate cases, Equation 24 becomes: 
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For low rate cases, Equation 28 becomes: 
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For very slow roll cases, Equation 30 becomes: 
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Comparing Cases 91,92, 93 in Table X with Cases 10, 11,12 in 
Table II and Cases 58,59,60 in Table VII it can be concluded that; 
except for possible computational error, Cases 91,92 and 93 are the 
algebraic sum of Cases 10,11, 12 and 58,59, 60; that is, for example, Case 
91 equals Case 10 plus Case 58. This fact is obviously true of the theory 
equations and is here shown to be the case for the 6-ID computations as 
well. Similar comparisons can be made with corresponding cases in 
Tables Ul, VIII, XI and H, IX, XU. Thus, the effects of aerodynamic 
asymmetries and those of initial transverse velocity are independent of 
one another, 


Figures 12,13 and 14 illustrate the Cases 91-123, The curves are 


me 


of the same form as Figures 7,8 and 9 but differ with th. addition of S,. 


Maximum effect of al) parameters is,desired. Cases 113, 114 and 115 show 


the limit of parameter combinations by 113 and 114 going unstable. Roll 
( rate effects are again large and velocity effects are larger than in Cases 


58-90. Figure 15 shows this to be true ii also shows the cases involving 
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TABLE X 


THEORY VALIDATION, ASYMMETREES, 
CASES 91-101 


Coefficients 


iA. (mils) 


m Eee on 
20.0114 20. 018+ 
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42, 325i] 43. 210i 
-9, 908 -}-10. 372 
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-29, 743i - 30. 427i 
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_ TABLE XI 


THEORY VALIDATION, ASYMMETRIES, 
CASES 102-112 
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TABLE XII 


THEORY VALIDATION, ASYMMETRIES, 
CASES 113-123 


Coefficients 


yA. (mils) 
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100, 3514] 100.0234 

115 6283 | 1100. 814i | 99, 986i 
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116 900 100. 587i| 99.7714 
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118 me QO} 0 100 } 1000 Al Al Al ieee cue wasi| 

102, 6684] 102. 349 
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i 
U = 3000 ft/sec to be ones or smallest dispersion. Such was the case in 
Figure 16. Figure 16 illustrates a sample trajectory, Case 101. 
Cases 124-156 
id 
To establish the relationship between the effects on dispersion for 
aerodynamic asymrnetries and initial angle of attack, a third set of cases 
were run, Again roli rate and velocity were varied as done previously but 
a, was set at (1+i) degrees with § =O and a. = 0. ‘Tables XIIf, XIV, and 
XV tabulate the results. For ali high roll rate cases, Equation 24 reduces 
to: 
a ip, oa pwd? ly- Ty i 
LS Aw 3(4 ~\y st a aor 1 
J.A.=) ag “80+ in ais M, : p “2,8 3, mud “+ “Fy |} 1000 
ie i. A 
For jow roll rato cases, Equation 28 reduces to: 
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rmx oie. € te M yy Ms‘ ¢ p" ( 
ae oe oa Sap 
ee - -= sin =—)}]1 
p \ U p sift «| CGO 
For very low roll retes, “quation 30 zuduces to: 
y 23 ™ 4 
A = PEOkt a [Cee ae | et x 2 i px \ 
i Deie e cae Weve A (lr Gore eo Gr 
L \ a e JL é 
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, f DX raG ry 4 (1% ‘, 
+ if $e sean \-— 
30 "4 29.20 aw 1000 
-Gnly for high roll rates does the &, term appear. @%, should have no 
noticeable effect on dispersion for p< 100 rad /sec. 
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TABLE XII 


THEORY VALIDATION, ASYMMETRIES, 
CASES 124-134 


— 
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TABLE XIV 


THEORY VALIDATION, ASYMMETKAES, 
CASES 135-145 
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i TABLE XV 


THEORY VALIDATION, ASYMMETRIES, 
CASES 146-156 


p 


31416 
18850 


O 
590 
300 
1. 

50 


25 


10 
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Coefficients 


J. A. Gnils) 


6-D | Theory 


Unstable 


Unstable 


0.0464} 0.008+ 
0.0391} O.00L1i 
0.237- 
0. 228i 
0.393- 
0. 384i 
1.2)3-}| 1.174- 

2.381-} 2, 349- 

fl 2.2941] 2.352% 

4.7\19-) 4,699~ 
4.637i| 4. 702i 


10,258-| 10.177- | 
14.4111] 14.4763 


24,440-| 24, 510- 
30.9761] 31.212% | 


-58.669 |[-58.450 
“58. 684i]-58. 450i 
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Comparing Cases 124,125,126 in Table XIII with Cases 19, 20,21 
in Table III and Cases 58,59,60 in Table VII it can be concluded that 
Cases 124,125 and 126 are the algebraic sum of Cases 19,20,21 and 58, 
59,60; that is, for example, Case 124 equals Case 19 plus Case 58. This 
is obvious from the reduced theoretical equations for Cases 124-156. It 
is shown here to be also true for the 6-D computations; allowing for some 
computational error. Similar comparisons can be made with correspond- 
ing cases in Tables III, VII, XIV and If, IX, XV. Thus the effects of 
aerodynamic asymmetries and those of initial angle of attack are inde- 
pendent of one another. 

Figures 17, 18 and 19 illustrate Cases 124-156. The curves are 
very similar to those in Figures 7,8, and 9 with the only difference being 
the very small a, contribution in Figures 17, 18, and 19. Cases 135, 
146, and 147 result in instabilities, indicating that maximum effect of the 
various parameters has been accomplished. Effects of roll rate are 
essentially the same as in Case 58-90 and effects of velocity, Figure 20, 
the same as in Figure 10. Cascs with U =3000 tt/sec again have the 


smallest dispersion. Figure 21 shows a typical trajectory, Case 134. 


Cases 157-189 
To validate the relationship between the cffects on dispersion for 
aerodynamic asymmetries and those of initial angular rate, a fourth set 
of cases were run. As before, roll rate and velocity were varied, but 


—b a — 
a, set at (250+ 250i) rad/sec with §,= 0 and a,=0. Tables XVI,XVII, 


184 


oe ae re ay 


| Y-FEET 
-60 -40 -20 0 20 40 
3 os 
t 
P=10 
oO uf 
: ° PaS 
| P=0 
a 
a 
ul 
uJ 
a P=2s 
m™ 
| o 
i is 
P=S0 
} 
| 
| | ras P=100 
| , Ss MH 
 e P=200 
ae | —  THEGAY P=S0Q 
ae | m 6-0 | 4-p=31416 AAD/SEC 
i 
| ol 
| ' N 
| Figure 17. Dispersion: Phase ILI Cases 124-134 
4 
1 4 ; 
Q : \ ‘ 
| : 185 
i : ‘4 
| ' 4 4 
y 7 ae ah ee A Dees au r 
an iraniielarotiea Aetairtiaciet WOH Seite nee oe pokmaskentenee eS alin eet ee ye 4 


med ae ee 
a cee Mae et net cet ak ok etn 


ae 
hee OT 
| 
| 
i 
/ 
| 
1 
| 
C. 
Cae : 
i 
; Meat 
' 
: 


. eee ee 
fs nares 


~60 


-26 


P=100 
P=i88S0O AAD/SEC 


aa THEQAY 
oO 6-0 


Figure 18. Dispersion: Phase III Cases 135-145 


186 


Lay 
arene waa (4 tenn 


4 
| 
i he 
a 
~60 
i) 
Po) 
’ 
Q 
@ 
’ 
I 
> 
re 1 
Lay 
lug 
ti 
8 
\ Mm 
L=) 
~ 
. ) 
| 
; =) 
f THEOAY 
6-0) 
oQ 
' ee | 
’ 
: | 
| 
' Figure 19. Dispersion: Phase III Cases 146-156 
‘ 1 
‘ | 
ae 


: 187 | 


Si ENTS e ree: 


a Ta ‘ 
1 
"| : 
| { 
a 
( Py ae ! 
— “> | 
bee. 2 ~69 
© 
: 1 Usi000 FT/SEC 
le 
w 
& : | 
URS000 FT/SEC 


~UG 


U=3000 FT/SEC 


Z-FEET 


-20 


Figure 20. Dispersion: Phase III Theory, Cases 124-156 


188 


H 


pov oe 


Race Oour- Zh 


ree Oouwr Ja 


* 1000 
1000 


oO o 
=) oO 
a ‘@ < 
ioe) 
et 
oe 
n 
oy 
G 
ro) =) a 
oO a B 
w 
| aad tr = 
WwW uJ = 
MW 
i 0 
<o xO . 
Oo 2 = 
= > be 
et) 
uM 
} 
5p 
ed 
ny 
°° 4 
oO (ae ] 
cq ba i 


ee. ee 


XVHI gives the results. For high roll rates, the governing equation 


becomes: 


Lewy = pumd™ —_ i ly — 
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For iow rol] rates, the governing equation: 
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For very slow roll, the governing equation: 
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Comparing Cases 157, 158,159 in Table XVI with Cases 28,29, 30 in 
Table ITV and Cases 58,59, 60 in Table VII, it can be concluded that Cases 
157,158, and 159 are the algebraic sum of Cases 28,29, 30 and 58,59, 60; 
that is, for example, Case 157 equals Case 28 plus Case 58. This is 
obvious .rom the reduced theorctical equations for Cases 157-189. Here it 
is shown to be true for 6-D computations also, Any discrepancy can be 


attributed to computational error. Similar comparisons can be made with 


corresponding cases in Table IV, VII], and XVII. Thus the effects of aero- 
iynamic asymmetries and those of initial angular rate are independent of 
one another 
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TABLE XVI 


THEORY VALIDATION, ASYMMETRIE 
CASES 157-167 


191 


S, 


6.576-| 6.707- 
BLE 14, 5624 
| 18.366-| 19.077- 

27.8271! 29, 000i 


-9.690-[- 10.283 
387 |-65.283i | 


.730i |-64. 426i 


.094 [-52,500 
O94} -52, 500% 


‘ TABLE XVII 


THEORY VALIDATION, ASYMMETRIES, 
CASES 168-178 
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TABLE XVII 


THEORY VALIDATION, ASYMMETRIES, 
CASES 179-189 


Theory 


Unstable 


Unstabie 


Unstable 
-5.015 | -5. 
~5.503i] -5. 709i 
“4, 884 
-5.572i 


-3. 191 
-7.892i 


l 
33. 158i) 36, 752i 


-61.894 |-63.990 
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Figures 22, 23 and 24 illustrate Cases 157-189. The curves are 
a 


similar to those in Figures 7, 8 and 9 but are displaced by the 5 


contribution, Cases 168, 179, 180 and 181 indicate that maximum effects 
of the various parameters has been achieved in ather stable cases. The 
effects of roll rate and velocity follow the same trends as those in Cases 
58-0. Figure 25 shows the effects of velocity for Cases 157-189, Cases 
with U = 3000 ft/sec exhibit the smallest dispersion. A sample trajectory, 


Case 189, is shown in Figure 26. 


Comparison: High, Low, Very Slow 


Roll Rate Theories 


In Cases 58-189 the High, Low, and Very Slow Roll Rate Pheories 
are validated for various initial conditions and parameters, The theories 
have been applied for certain ranges in roll rate and roll rate times time 
(pt.) The range of pt, (pt <1.0) are governed by the inherent requirements 
of power scries expansion. However, the ranges of p are arbitrary (to a 
certain extent) and are based on accuracy of the theories themseives. Each 
theory approximates the solution very well for a certain range of p and then 
begins to diverge and become inaccurate. ‘The range of p for which the 
very slow roll rate theory is accurate is fairly well cut and dried; p20, 
pt 31.0, For any pt> i.Q we must now use the low roll rate theory, The 
question now arises, how high a roll rate can this theory accommodate? 
At what value of p must we change to the high roll rate theory? These 


questions are answered by a plot of sample 6-D computations, Figure 27 
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and all three theories éenaed beyond the limits used in the previous 
validation. The high roll rate theory is a straight line going off to infinity 
as p goes to zero. Although the length of the curve in which it is an effective 
theory is short graphically, the range of roll rates it encompasses is 
tremendous, Figure 28 illustrates the effective limits of each theory; that 
is, onthe spectrum of possible roll rates it shows where each theory is 
the most effective. The low roll rate theory handles the largest graphical 
area but only roll rates less than 100 rad/sec and greater than 5 rad/sec. 
The upper limit of 100 rad/sec was chosen since here the low roll theory 
attaches itself to the 6-D results while the high roll theory diverges. The 
lower limit of 5 rad/sec corresponds to pt<l.0. Figures 27 and 28 depict 
Cases 58-68 where uy = 5000 ft/sec or t = 0.2 sec, Therefore p = 5 rad/ 
sec corresponds to pt=1.0. The very low roll rate theory has the smallest 
range but is essential in predicting dispersion as the roll rate goes to zero. 
As pt>lL, the theory diverges as would be expected from a power series; 
Equation 29, The sharp turn occurs at p&20 rad/sec or ptx4 for Cases 58- 
68. Although Cases 58-68 were illustrated here, this analysis of the 
effective limits of the roll theories was found to be similar for all other 
cases, For the u)= 3000 ft/sec cases the iow roll theory limits were 


3,0<p <S0 for tee 1006 ft/sec cases: 1.0<p<25.0, 


Phase IV 
To validate the effects of gravity on dispersion, a final set of 


cases were run using the high roll rate theory, Equation 24, Ordinarly, 


one would think that gravity would only introduce a constant term; one 
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Figure 27. Phase III Theory Equations 24, 28, 30 
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that could be factored out. However intcgration of the equations of motion 
produce a eracity term dependent upon roll rate. Deterrnination of its 

validity and consequence is what is important here. S <a, and @, were 
set to zero in crder to allow determination of the effects due to roll rate 


and velocity. The reduced governing equation becomes: 


JA.= ig =r) + ws 1000 


No acrodynamic asymmetrics were present and the effects of gravity were 
assumed independent of effects due to 3, Qo; a; a logical assumption, 
Table XIX lists the results. 

Table XIX indicates that the effects due to gravity occur largely 
in the vertical plane, as would be expected. The transverse contribution 
is minimal but is affected by both velocity and roll rate. The vertical 
contribution is only affected by velocity. The unstable cases indicate 
maximum use of Magnus and thus maximum transverse effects on dispersion. 
It can be concluded from this brief but thorough treatment that gravity 
effects dispersion only in the vertical plane (for all practical purposes) 
and that its contribution is constant with velocity. The roll dependent 
term, ipkk , Has bee shown to exist but becom: negligible for the 

mud 

flechette. This term would possibly become important for projectile dis- 
persion and other missile applications. Projectile motion with gravity is 
typified by a cocking right of the proicctile in flight with a positive OM y 


but negative Cy? the parameter A would become negative and the entire 


roll dependent term, positive; that is, cocked to the right, dispersion to 
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TABLE XIX 


THEORY VALIDATION,GRAVITY 
CASES 190-201 


C Initial Conditions | 
| a | 
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the right. For a finned missile the opposite would occur due to the 
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FREE FLIGhT DATA ANALYSIS 


In order to analyze actual test firings as to jump and dispersion and 
correlate them with the validated theory, the initial conditions of each 
test firing must be obtained and ptit into the proper form. To obtain raw 
experimental data, test firings were conducted by the U.S. Army, 
Frankford Arsenal. The configuration tested was the Producibility Ground 
Point Flechette, Figure 29, The raw data required was both translational 
and angular; that is, data was needed to determine position as a function 
of time and angle of attack of the flechette as a function of time. To 
accomplish this, Frankford Arsenal devised the test apparatus shown in 
Figure 30. The gun barrel was mounted on a steel girder and a laser 
beam was used to obtain the aim point on a target 50 meters down range. 
At positions, 1, 3, 5, 7, 9, and 11 feet downrange, orthogonal flash x- 
ray tubes were placed to photograph the flechette as it passed its station. 
One tube was placed to allow a top view at each station and provide a 
means of obtaining swerve and yaw data, The other tube allowed a side 
view at each station to obtain heave and pitch data, At each station reference 
marks oriented the flechette as to its exact position downrange. This was 
to allow for any timing error and/or variation in muzzle velocity. The 
photographs were taken using special soft flash x-ray tubes which permit 
the photographing of the low density sabot pieces and analyzing the separa~ 


tion in addition to the mction of the flechette. 
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From the battery of test firings of 20 rounds of each type (which 
included tests of the ground point, and swayed point producibility flechette 
as well as the R&D version), 8 of the ground point producibility rounds 
were selected to be analyzed. The eight rounds along with velocity, roll 
rates and target positions are given in Table XX. 

Raw transiational and angular data are shown in Figures 31 through 


46. The figures illustrate the position and complex angle of attack of the 


flechette for each station, 


TABLE XX 


FRANK FORD TEST FIRING DATA 


Target at 50 ft. 


p Y (ft) iZ. (ft) 


(@) 
11,454 
13,201 


14,219 


053 | 0.099 
13,289 0.053 0.016 
17,304 | 0.084 -0.004 | 
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Figure 35. Raw Translational Data Ground Point - Round 7 i 
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Figure 36. Raw Angular Data Ground Point -— Round 7 
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Once the raw data was obtained, it had to be converted into a form 
such that initial conditions So» So » @ and a, could be extracted from it. 
as ce 
To eventually arrive at values for So and So , the translational parameters, 


the raw position or translational data had to be approximated by eqinations. 


The raw data was fitted to a polynomial equation of third degree by a least 


squares method. The data in the y-direction was fit separately from that 
in the z-direction to distinguish between the swerve and heave contributions. 


With the equations obtained, a simple differentiation yielded equations for 


od 
i e ‘ 
the velocities in the y and z directions. The initial conditions So and So i 


are now rcadily obtainable: 


—b 
35 (ft) = ye + iz, 
aad 
So (it/sec) = Yo iz, 
é = eres 

Obtaining a, and &, from the raw angular data was more difficult, 
The traditional way of analyzing any missile motion with pitch, yaw, 
and roll is by a three-degrec-of-freedom lcast squares fit to the tricyclic 
motion, Equation 6. Hlowever, the availability of only 6 data points made 
this technique impossible, so another, approximate mcthod, had to be 
employed. The solution was to approximate the pitching and yawing motion 
to one -degree-of-freedom while holding the roll rate constant. In order to 
do this, the B - @ axis system had to be rotated to coincide with the more 
dominant angular mode. Figure 47 illustrates a typical raw angular data 


plot. Since the angular motion of the flechette tends to approximate an 


cllipse, the B -a@ axes are rotated some angle y to coincide with the 
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Figure 47, Axis Rotation Approximates Pure Pitching Motion ig 


| : 


major and minor axes of the ellipse, 4s shown. The angular data is 
retabulated for this new axes system, p'- a’. To fit the data to the one- 
degree -of-freedom equation: 


a =Ke*t cos (wt +6) 


only the dominant mode can be considered, For example, in Figure 47 the 
dominant mode occurs along the a" axis; therefore, only @' coordinates 
are utilized in the least squares fit, corresponding g' coordinates are 
ignored, Table XXI lists the parameters obtained for the eight flechette 
rounds. Once an equation for a@' is obtained, it represents one dimensional 
oscillatory motion along the q' axis. A simple differentiating of the a‘ 

She 


equation yields an equation for a@', The initial conditions @ 
O 


=e 
and @,, 
however, are complex whereas q@' and @' are only one dimensional. 
Therefore, the rotation angle y is taken into account and the a' equation 


is projected back into the B, @ axes system: 


i 


a=a' cosy 


a’ cosy 


i} 


Bp =a’ siny 


Tous the complex initial conditions are approximated. 


@ = Bp. + ia 
% = Bo ) 
~ ‘ | | 

. 2? { 
a, = 6, + id, | 


Figures 48-63 illustrate the fitted data both translational and anguiar 


for the eight rounds. The transitional data includes the pertinent equations. 
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Figure 56. Fitted TrausJuttonal Data Ground Point - Round 14 | 
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Figure 57. Fitted Angular Data Ground Point — Round 14 
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Figure 58. Fitted Translational Data Ground Point — Round 16 
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Figure 61, Fitted Angular Data Ground Point - Round 17 
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Figure 62. Fitted Translationai Data Ground Point - Round 19 
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Figure 63. Fitted Angular Data Ground Point - Round 19 
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DISPERSION ANALYSIS 


Free Flight vs Theory 


Once the initial conditions are determined as in the previous section, 
they are applied to the theory and compared to the dispersion of each test 
fired round. To utilize the theory, the fitted data must be chosen for a 

oak Sy ae 
given time; that is, So: So Us % Must be selected for one given point 
in time - position dowirange. Since the question of what poirt in time do 
the initial conditions occur, 3 sets of initial sGHaiions were chosen to 
correspond with positions 1, 3, 5 fect downrange. Vhis span of position 
downrange may or may not be sufficient to include the actual time correspond- 
ing to the initial conditions for 2ach round. The following analysis will 
determine each round's effective time for its initial conditions. 

Por each set of initial conditions, theory and 6-D computations were 
done and compared to target data for the Frankford test firings. The results 
are tabulated in Table XXII in mils and plotted in Figures 64-71 in feet; 
deviation from the time of fire at 50 ft. downrange. The relationship 
between the deviations in feet and mils at 50 ft, downrange is: 


=e 


J.-A, (mils) = S (ft) (1000) 
x 


or J-A, (mils) = (20) 'S (ft) 

Yo accurately and concisely analyze the complex and large amount 
of data in Table XXII, the positions downrange in which the initial conditions 
were selected must be simultaneously analyzed with the aispersion results 
at 50 ft downrange. The problem in choosing initial conditions is where 


they should be taken; at what point downrange. Normally, one would think 
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ae 


that the initial conditions would occur immediately after leaving the gun 
barrel. However, the flechette being a finned body needs a sabot con- 
figuration to guide it down the barrel, Figure 29, The sabot causes the 
initial condition location problem since the sabot must separate from the 
flechette outside of the gun barrel, The exact time and place where this 
occurs is not constant; varying from round to round. Not only does the 
sabot separate from the flechette instantaneously different every time, 

the sabot may not separate cleanly or the same way every time. Inter- 
ference with the fins after sabot separation can cause disturbances to the 
flechette and alter the initial conditions, In addition, asymmetric sabot 
separation can influence the initial conditions. [Tigures 72-79 illustrate 
the flight transition sequence for the 8 flechette test rounds, In every 
sequence the sabot begins to separate, in varying degrces, | ft. downrange. 
At 3 ft. downrange,the sabot is nearly completely separated, but in some 
cases the sabot particles pose interference problems with the fins. By 5 and 
7 ft. downrange the sabot has completely separated and the flechette is in 
free flight. The correspondence between the flight ecnction sequence and 
dispersion results can be seen in each individual round. Figure 64 indicates 
that the initial conditions for round 4 occur somnewhere between | and 3 ft. 
downrange judging by the dispersion of the actual tested round. Figure 72 
verifies this fact in that the sabot has separated from the flechette between 


i and 3 ft. dcwnrange. The y-coordinate in the dispersion vector does not 
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Figure 79. Flight Transition Sequence - Round 19 
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accurately agree with the theory for this case, However, besides 
| computational error other physical factors can influence dispersion. 


Contributions by fin asymmetries and other configurational asymmetries 


can be important but are unable to be detected or accounted for. Through - 
out this analysis this must be kept in mind to partially account for any 
discrepancy between the actual test firing and the theory and 6-D 
computations. Figure 65 indicates the initial conditions for round 6 occur 
between 3 and 5 ft. downrange. Figure 73 verities this choice showing 
separation occurring around 3 ft but with sabot particles very close to the 
Ls fins causing possible interference and delaying the initial conditions 
location. The initial conditions location for round 7 is difficult to accurately 
choose since the y-coordinate does not accurately agree, igure 66. It is 
safe to say that the initial conditions occur sometime around 3 ft and 
Figure 74 vcrifics this choice. The %-coordinate for round 8 is not as 
accurate as would be desired, Figure 7, but the y-coordinate indicates 
initia! conditions occurring between 3 and 5 ft downrange. Figure 75 agrees 
with this choice indicated interference with the fins at 3 ft deiaying the 
initial conditions. Initial conditions for round 14 are chosen between 3 and 
5 ft. downrange, igure 68. Figure 76 indicates possible fin interference 
tending to verify the choice. Figures 69 and 77 indicate and verify the 


i 
| 
Y 
t 
{ 
|: 
‘ choice of initial conditions in the inimediate vicinity of 3 ft downrange for 
| round 16. Possible fin interference at 3 tt downrange, Figure 78, round 17, 
| 
t 


; verifies a choice of initial conditions between 3 and 5 it, Figure 70. A similar 
| : ds situation occurs for round 19 in Figures 71 and 79. Jt is often difficult to 
! y 
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choose initial condition positions accurately due to slight discrepancies 
between theory and test firings. However, the discrepancies are of the 
order 0.05 ft, which shows up large in Figures 64-71 due to the scale 
chosen, but is within the error expected from the validation of theory 
section. 

The influence of sabot separation can be readily seen by inspection 
of Figures 72-79, 1 and 3 ft downrange. In every case, the flechette and 
gabot are at nearly a zero angle of attack at 1 ft, but has changed angle of 
attack noticeably by 3 ft downrange. This would indicate that fin inter- 
ference or asymmetric sabot separation is causing the noticeable effect. 
It can be concluded that dispersion is dependent upon the initial conditions 
that the initial conditions are a function of sabot separation and that the 


theory can predict what the initia] conditions are and where they occur. 


Dispersion Theory vs, First Maximum Yaw [Hypothesis 

A popular theory to predict the dispersion of flechettes is the First 
Maximum Yaw Hypothesis. This theory relates the digpersion magnitude 
to the first maximum yaw magnitude by a nearly linear relationship. Other 

——_ ~ 

initia! conditions such as angular rate, @, and transiationai velocity, Sg 
are said not to effect dispersion. To disprove this theory and strengthen 
the position of the theory ascribing to dispersion due to initial conditions 
® ee “a ar 
So: ayy b> the First Maximum Yaw theory was applied to Frankford 


Arsenal data, Figure 80 shows a plot of dispersion magnitude vs. first 


maximum yaw magnitude. Clearly no linear relationship exists between 
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dispersion and first maximum yaw. In fact, the plotted data resembles a 

| random shotgun blast, Figures 81,82 and 83 employ the theory to the first 
maximum yaw hypothesis. Again the plot substantiates the findings of 

Figure 80. The disproval] of the first maximum yaw hypothesis comes as no 
surprise since the dispersion theory contradicts it and the 6-D computations, 
which integrate the actual equations of motion, validated the dispersion 
theory. Therefore, dispersion could never accurately be predicted by a 
theory involving only first maxirnum yaw. 


obs 


Ld — « 2 
The influence of initial conditions, So» a , and a, and dispersion 


{> for these rounds. 


for the actual test firings are expected to be different from that in the 
. validation of theory section because of the different ranges in the initial 
: — 
conditions. For example, 5, in the validation section was (100 + 100i) 
=k 
| ft/sec. In the actual test firings, Sy only ranged a3 high as 0.017 ft/sec, 
ed 

| Of course, the large value was only to validate the theory. Here Sy is very 
| ( small and its contribution is accordingly smaller. In the reduced equation 
“| 
a 24, employed to calculate the theory colunin in Table XXII, 

—_ => 

: ils -=l 20 . Ras - — - x 
3 J. A. (mils)= 1000 : 4 7 da (&, ay hy 

/ 
roi for round 4, | fe downrange, 
. : 
| ‘ 1000 —" = (0.001562 + 0.001841i) mils 
| : where as, 
| _ 
| | J.A. = (1,329 - 1.302%) mils 

| 3 ra 
' Since this is typical of the 8 rounds tested, 5, has little effect on dispersion 
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Similarly, for this particular case, 7 


oven 


S 
1000 — = (1.519360 + 0.041760i) mils 


| 
ipL.A 
1000 a JPA. (.0,01437 + 0.002141) mils ? 
Oo mud 
e Obviously, SS and a, are by far the greatest contributors tc dispersion 


for this case. Inspection of all the other 23 cases in Table XXTI agrees 


ae. 


R with this general pattern. So can be nearly eliminated, of course, by 
accurate setup of the test equipment so that the gun barrel is set exactly 
at coordinates (0,0). Any 5, then would occur from displacement due 
to the blast. This leaves the major culprit in dispersion to be & ; 
Figure 84 illustrates the dependence of the Jump Angle, and hence 
dispersion, upon angular rate and angle of attack. 

Although @, contributes the most to the Jump Angle, the combination 


te 


of So and a also has a noticeable influence. From the test firings, S, 


gor 


was found to have a negligible effect on dispersion. Therefore, it is 


neglected in Figure 84 to simplify the plot. It is evident from Figure 84 


S eee ee ae 


—ae 
2 = - oly, e : 2 
that various combinations of a, and a, yield zero dispersion. It is 


cn 


possible that large values of a and &, can combine to yield zero dispersion; 
an impossibility with the first maximum yaw hypothesis. If ay and @, are 
able to balance to give zero dispersion, then this idea can be expanded to 


f include the entire equation. 
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The governing equation used throughout this dispersion analysis 


section is: 


—_ ele 


cae So So I 
J.4.= 1000) = + - + a(& “Ge =) + ( 7) 


Eliminating the constant gravity term, 


Setting J. T.A. to zero, the idea benind Figure 84 is expanded to include 


te * 
So» So- 
oD So = A & zo Phx 
x u mud ‘0 % 
rearranging 


— u aris A — Ses 2 
m SoH + Sol= =p (4% Yo ipl.) 
A dimensional analysis of the equations finds that both sides have units 


of momentum or impulse. Going one step farther it can be said that to 


obtain zero dispersion: 


initial transverse momentum-= initial angular momentum 


Therefore it is the imbalance in the initial mo 


inentums that causes 
dispersion, The size of initial conditions can be huge, Figure 84, but 
if they can combine to balance, zero dispersion results, The way the 
initial conditions combine, determine the magnitude of the imbatance or 


dispersion. It should be noted that this dispersion discussed is round to 


round dispersion and that the inconsistency of the momentum imbalance 
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from round to round causes a dispersion pattern (a set of rounds). The 
next section will highlight this principle in the evaluation of physical 


factors affecting dispersion. 
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PHYSICAL, EVALUATICN OF DISPERSION 


Initial momentum imbalance has been shown to cause dispersion. 
Initial conditions determine the magnitude of the imbalance. What causes 
these initial conditions to occur is the subject of this final section. Initial 
conditions occur somewhere between zero and five feet downrange to 
different degrees of magnitude due to various conditions. These conditions 
are: 

1. Fin or body asyAieny 
2. In-bore mal-alignment 
3. Asymmetric blast 

4, Asymmetric sabot separation 

5. Sabot-fin interference 

6, Fin or body damage 
Fin or body asymmetries can cause dispersion magnitudes to range as 
much or greater than those in the Validation of Theory section for aero- 
dynainic asymmetries. These asymmetries can be avercanted or bent 
fins, damaged nose cone, or even body deformities. Figure 85 which 
shows in-bore imal-aligninent aiso shows a slightiy bent body, concave 
downward. In-bore mal-alignmert can be attributed to warping and/or 
the entire flechette at some angle of attack. Clearly, if this flechette were 
fired, the in-bore angle of attack wouid produce an aX, outside the gun 
barrel even before sabot separation, With the flechctte at some angle of 


ee 


— —» ° 
attack, the blast can cause a large dy and an So and So . the blast itself 
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is a chief catalyst in causing the initial conditions. An asymmetric blast 
can indeed impart influence on the initia! conditions, bt ymmetric 
blast can also. Given an initial angle of attack due to sor tisturbances 
the symmetric blast can cause significant a, > 8 and 8. Figure 86 
shows a typical blast region with the flechette outlined in the picture. 
The momentum principle discussed in the previous section goes hand-in- 
hand with this blast region. It is here that the transverse and angular- 
momentum is imparted to the flechette. Figure 87 illustrates a typical 
flechette in the blast region. Coming out of the barrel at some angle of 
attack, the blast catches the flechette and induces some angular rate. 
At the same time, the flechette is translated laterally giving an Sy and 
s If these contributions cancel each other out; that is, if initial 
transverse momentum equals initial angular momentum then the dispersion 
is zero, If they do not cancel, dispersion results. The sketch is highly 
simplified in that the blast itself is all-engulfing as in Figure 86, Of 
course, the transition sequence of sabot separation, fin interference, 
and possible fin damage must not be forgotten. The transition sequence 
occurs in the blast region, however, and is not considered separate 
from the blast. When separation occurs, the sabot particles are apt to 
interfere with the fin section and cause possibie damage. Once the sabot 
has separated and cleared the fins the blast has had its greatest effect 
and the initial conditions can be determined. After the flechette has 


moved downrange, it assumes supersonic free flight, Figure 88. 
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Free Flight, Ground Point Flechette 


Supersonic 


Figure 88. 
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CONCLUSIONS 


A complete Jump and Dispersion Theory has been developed for 
free flight vehicles. Three governing equations have been determined to 
accommodate high, low, and very low roll rates. The theories were found 
to be accurate with six-degree-of-freedom numerical computations of the 
equations of motion and therefore reliably predict the jump and dispersion 
of flechettes. The theory validation included 201 case runs in four phases. 
The first phase validated the theory with respect to restoring and damping 
moments. The effect of these moments on dispersion was found to depend 
on the initial conditions. The second phase validated the theory with respect 
to Magnus forces and moments. The effect of Magnus was found to be very 
small and not to be of any consequence unless the total dispersion of any 
given round was of the same order of magnitude as the Magnus effect. 
Phase three validates the theory with respect to aerodynamic asymmetries 
and roll rate. All three theories were validated in this phase and found to 
be quite accurate considering the large dispersions encountered. Aero- 
dynamic asymmictries Causing a trim angle of 1° had little effect on the 
dispersion of flechettes. Slower rolling bodies were shown to have, in 
general, increasingly larger dispersion values as roll rate aeereaeey. 
It can be concluded that for free flight vehicles that are prone to aero- 
dynamic asymmetries and fin damage, a high roll rate is essential to 
lower dispersion and increase accuracy. The fourth phase validates the 
theory with respect to gravity. The theory indicates a lateral contribution 


to dispersion from gravity in addition to the obvious vertical contribution. 
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For the flechette, the lateral contribution was found to be minimal and 
was neglected in this analysis. 

Free flight data was obtained from Frankford Arsenal to correlate 
with the theory. Angular and translational data was fitted and put into 
initial condition form. The initial condition data was applied to the theory 
and compared to target data for the rounds tested. The theory was found 
to agree favorably in magnitude with the test firings. As a result, the 
method used to analyze the data can be considered a valid method. Photo- 
graphs of the test firings were taken to include the flight transition 
sequence in the blast region. The pictures further verify the analysis 
method of the initial conditions by allowing agreement between the chosen 
initial conditions and the position downrange where they were selected. 

The evaluation of the free flight dispersion against the theory also 
disproves the First Maximum Yaw hypothesis. A plot of jump angle vs. 
first maximum yaw of actual test data produced a shotgun biast pattern 
with no relationship evident between dispersion and first maximum yaw. 


. 


In addition, a plot of jump angle versus angular rate for various initial 
angles of attack indicates an infinite amount of combinations o 
conditions to yield a given jump angle. Thus, zero dispersion has an 
infinite set of possible initial conditions. It was found for zero dispersion 
that a unique physical condition holds: to obtain zero dispersion, initial 
transverse momentum = initial angular momentum. These impulses are 


imparted to the flechette in the blast region where the body and especially 


the fins are subject to disturbances. Momentum imbalance is the reason 
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dispersion occurs, The initial conditions only determine the magnitude 


of imbalance or dispersion. This dispersion is round to round dispersion. 


Inconsistency in the imbalance results in a dispersion pattern. The 

initial conditions were found not to occur until after the sabot separation 
and the blast has had its greatest effect. The factors causing the existence 
of initial conditions were found to be not only the blast and sabot separation 
sequence, but also fin and body asymmetries and bore mal-alignment. 

In order to decrease dispersion, these physical factors causing initial 
conditions must be kept at a minimum. The most important aspect would 
be to protect the fins from asymmetries, damage, and interference from 


the separating sabot. Initial conditions can never realistically be 


eliminated but if kept minimal, dispersion is reduced. 
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Appendix Al contains mass parameters and stability coefficients 

for the Ground Point Flechette. Table Al-1 lists values for mass, diameter, 

| axial and transverse moments of inertia. Figure: Al-1 through Al-8 
present stability coefficients used in this analysis versus Mach number. 
Czy CMY? CMg + CM, were provided by Frankford Arsenal, Czpp ’ 

My, » Cyp, Czp, CMg» Cryg were nominal values of the coefficients 


following the same trends of Czy and CM for Mach number. CM g and 


. | CM, CM 4 were verified in the University of Notre Dame supersonic 


wind tunnel. 16 
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TABLE Al-1 


FLECHETTE PARAMETERS 


mass 
diameter 


L, 


I 
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0.000046 slugs 

0.006 ft. 

0.000000000217 slugs-it? 
0.000900036421 slugs-tr* 
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Figure Al~l, CZa vs Mach Number Producibility Ground Point 
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Figure Al-8. 
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APPENDIX 


A-2 


Appendix A-2 contains the complete print-out of the results from a 
typical 6-D computer program run. The results give the time from launch, 
position coordinates x,y,z, velocity, roll rate, the magnitude of the 
complex angle of attack, Mach number, roll orientation angle, angles of 
pitch and yaw, nutation and precession damping factors, nutation and 
precession mode frequency rates, the gyroscopic stability factor, dynamic 
weight factor, and trim angle, 

The program is divided into various subroutines to eliminate any 
superfluous calculations, These subroutines read in aerodynamic ccefficients 
in tabular form as functions of Mach number and angle of attack, initialize 
the data, and integrate the six-degrees~-of-freedom differential equations 


of motion using a four-step Runge-Kutta scheme to obtain the vehicle 


trajectory. 
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FIRING PROGRAM OF FLECHETTES 


During the spring and summer of 1974, eleven flechette firings were 


og . carried out in the Frankford Arsenal X-nay Ballistic Range. The data 
‘ # a i 
=| was nleasured from the x-ray plates by Frankford Arsenal personnel and 
ae is given in Figure D-1,. 

ae 
| : The flechette used in this test is very similar to the flechette used 
~~ oe in lot 3 of the previous test data. It is assumed that the C.G. and the axial 
+ alt 

¥ 
ark and transverse moments are the same. The C.G. is 0.94 inches from the 
: é nose; the axial moment is 0.0107 grain in? and the transverse moment is 
wer . : 
* sie 1,78 grain in, It is also assumed that the average spin rate is 2750 
A eee 9 Te rev /sec. 
ny fe The coordinate axis for the target data is defined as having plot (0,0) 
4 a irs 

ae? coincident with the aim point as established with a laser. 

er Times of flight were obtained between the first x-ray station and the 


sixth x-ray station. The baseline for determining the velocity was then a 
variable where the distance was measured between the flechette C,G. at 
the first station and at the sixth station. 

The data used in the ND analysis is given in Figure D-2, The horizontal 
and vertical angles for Round 17 are plotted in D-3. In the case of Round 17 
and Rounds 35, 23, 21 and 25 the motion was nearly ld and therefore 


reductions and analysis could be carried out. The motions, however, on the 
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other rounds were highly 2d and therefore computer fits were not possible. 
The results for Round 17 are presented herein in order to provide the reader 
with an example. Figure D-4 contains the original a, g data and the new 

data as determined by computer fit. The finala, g data is given in Figure D-5 °° 
which is used in the basic jump equation. 

The trajectory data given in Figure D-1 is fitted on the computer 
with second degree, third degree and fourth degree polynomials as shown 
in Figure D-6. The third degree fit was used and a summary of the 
trajectory data is given in Figure D-7 which is also used in the basic jurnp 
equation. 

i Therefore, the @ andg data from Figure D-5 and the trajectory 
data from Figure D-7 is used with the basic jump equation to provide the 
ides results given in Figure D-8, 

Figure |--8, thercfore, is the final result of the firing program and 
the jump analysis which shows the agreement between the experimentally 
determined jump and the jump predicted from the translational ballistics 

’ \ 
motion parameters as determined in the Frankford Arsenal X-Ray Ballistics 
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R&D PARAMETERS 


d= 0.00587 ft. | 

m = 0Q.00004573 slugs 

L, = 0.000 000 000 330 slugs-tt? 
= 0.000 000 054 883 slugs-tt* 


p = 17279 rad/sec 
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